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Abstract. We study sublcading corrections to the corner free energy in classical two-dimensional critical 
systems, focusing on a generic boundary perturbation by the stress-tensor of the underlying conformal field 
theory (CFT). In the particular case of an angle 27T, we find that there is an unusual correction of the form 
L logL, where L is a typical length scale in the system. This correction also affects the one-point function 
of an operator near the corner. The prcfactor can be seen as semi-universal, in the sense that it depends on 
a single non-universal quantity, the extrapolation length. Once this ultraviolet cutoff is known, the term is 
entirely fixed by the geometry of the system, and the central charge of the CFT. Such a corner appears for 
example in the bipartite fidelity of a one-dimensional quantum system at criticality, which allows for several 
numerical checks in free fermions systems. We also present an exact result in the XX and Ising chains 
that confirms this analysis. Finally, we consider applications to the time evolution of the (logarithmic) 
Loschmidt echo and the entanglement entropy following a local quantum quench. Due to subtle issues in 
analytic continuation, wc find that the logarithmic term in imaginary time transforms into a time-dependent 
L~ 2 correction for the entanglement entropy, and a L~ 2 logL term for the Loschmidt echo. 
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1. Introduction 



1.1. Universality and free energies 

Phase transitions in statistical models at equilibrium are among the most studied phenomena in physics, with 
many important results that have accumulated over the XX t?l century. Our theoretical understanding of these 
phenomena largely relies on lattice models, some of which may be solved analytically [1] — or at least studied 
numerically — as well as on field theoretic results [2]. Critical points fall into different universality classes, 
which are caracterized by the large scale properties of the correlation functions of local observables. The 
free energy of these systems can also reveal interesting information about their underlying universality class 
[3, 4], or about the renormalization group flows between different critical points [5, 6, 7, 8]. Quantitatively, 
our understanding of critical phenomena is most advanced in two dimensions, both from the lattice and the 
field theoretic points of view, thanks to lattice integrability [1] and conformal invariance in the continuum 
[9]. 

An interesting development over the past twenty years has been the application of the techniques 
initially developed in the field of classical statistical models to quantum many-body systems. There, some 
quantities that are typically of interest are the ones that characterize the entanglement in the ground-state. 
Prominent examples of such quantities are the von Neumann entanglement entropy, the Renyi entropies 
[10, 11, 12], or other types of fidelities [13, 14, 15, 16]. In many situations, thanks to the classical-to-quantum 
correspondence, such quantities, defined for the ground state of a quantum system in d spatial dimensions, 
can be reformulated as free energies of a classical system in d + 1 spatial dimensions [10, 12, 15]. It turns 
out that these classical statistical systems are usually defined on domains that include conical singularities, 
and/or corners if the systems have boundaries. For the convenience of the reader, some of these aspects are 
reviewed in the Appendix A. 

Motivated by these recent developments, the purpose of this paper is to study classical statistical models 
with a boundary, in the presence of a corner. For critical two-dimensional systems described by a conformal 
field theory, it has long been known [17] that the free energy contains a piece of order 0(log L), where L is the 
typical size of the system, with a universal prefactor proportional to the central charge. This phenomenon 
has quantum-mechanical counterparts: in [15] we showed that the same universal term appears in the 
(logarithmic) bipartite fidelity (LBF in short) of the ground state of a quantum critical system with dynamic 
exponent z = 1, or in non-equilibrium situations, after a local quantum quench (see also the Appendix A). 
We will come back to quantum systems and bipartite entanglement later; for now, let us formulate our 
arguments in the language of classical statistical systems at equilibrium. 

The study of statistical systems in the presence of corner singularities is a topic that is interesting in 
its own right [17, 18, 19, 20, 21, 22, 23, 24]. In this paper, we focus on the subleading corrections to the 
corner free energies in two-dimensional critical systems, beyond the leading order of Cardy and Peschel [17]. 
In general these subleading corrections are generated by irrelevant bulk or boundary perturbations to the 
action of the conformal field theory (CFT), and they take the form of 0(L~ a ) terms, with some exponent 
a (not necessarily an integer). Here we focus on the arguably most generic boundary perturbation, which 
corresponds to the stress-tensor of the CFT itself. When the corner has an internal angle 2n, and only in 
that case, we show that the stress-tensor generates an interesting subleading contribution to the corner free 
energy of order 0(L~ l logL). A similar correction also appears in the correlation functions near the corner. 
Such corrections may be regarded as semi-universal: they are of the form 

(^corner = C / goom . - bg(X/a), (1) 

where c is the central charge, / geom . is a dimensionless geometric factor that is easily calculable (we will give 
a closed formula for this factor below). The only length scale that appears here, apart from the typical size of 
the system L and the UV cutoff a (lattice spacing), is the extrapolation length £. This extrapolation length is 
a concept that is ubiquitous in the study of surface critical phenomena [25] . The importance of this length has 
been emphasized lately in the context of global quantum quenches [26] , entanglement spectra in fractional 
quantum hall systems [27], or impurity problems such as the Kondo screening cloud [28]. Its appearance in 
logarithmic corrections to the corner free energy provides a practically reliable way of measuring this length 
in numerical simulations; it also affects various computable quantities in one-dimensional quantum systems, 
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such as the bipartite fidelity [15] or the emptiness formation probability [29]. Below, we will also explore 
the consequences of the appearance of this correction in time-dependent problems following a local quench, 
and we will provide extensive numerical and exact computations in support of our results in free fermion 
systems. 

1.2. Why is 0(L~ 1 \ogL) interesting? 

A simple way of seeing why such corrections are robust is to write down an expansion of the free energy, as a 
function of the dominant length scale L. In the following we will encounter several times the following form 

F(L) = f 2 L 2 + f x L + b log L + f Q + 6-i ^ + ^ + o(l/L). (2) 

f 2 and /i are the bulk and surface free energies in two dimensions. They are sensitive to the microscopic 
details of the lattice model. Universal properties, associated to the correlations at large scales, appear in 
the subleading corrections. Typically, in applications to entanglement measures in quantum systems, the 
contribution of the extensive part of the free energy, f 2 L 2 + fiL, disappears because of the normalization 
of the reduced density matrix [12, 15]. This is the case both for the entanglement entropy and for the LBF. 
Then, in these applications, the first term of interest is the one proportional to logL. This logarithmic 
term has a dimensionless prefactor bo which does not depend on any cutoff: it is universal. Of course, non- 
universal short length scales do appear when one computes the free energies of lattice models; their possible 
effect can be guessed by formally replacing the dominant length L by L + a in all the terms! in (2), and 
expanding again in L. We get 

F' (L) = f 2 L 2 + (/i + 2af 2 )L + b log L + f + af x + a 2 f 2 + 6_r ^ + + {l/L) . (3) 

We see that the coefficient bo remains unaffected, as expected. The same thing happens for the 0(L _1 logL) 
term. For this reason, even though &_i has the dimension of a length, it is a quantity that is more robust 
than /i, fo or More precisely, b-\ must be of the form 6_i = £ x (universal factor), for a single length 
scale £. Once £ is fixed, then the coefficient b-\ is uniquely determined by the global geometry of the domain. 
As mentioned earlier, we will see that £ is nothing but the extrapolation length [25], and the universal factor 
/geom. = 6_i/(£c) is calculable with CFT methods (here c is the central charge). Of course the argument 
presented here is qualitative; it needs to be supported by field-theoretic as well as lattice calculations. This 
is the goal of this paper. 

1.3. Organization of the manuscript 

The paper is organized as follows. In Sec. 2 we explain the generic appearance of a L _1 logL term when 
there is a corner with angle 2tt. We provide explicit and general formulae for the free energy as well as 
the one-point function of a primary operator. We then apply these result to the "pants" geometry that 
corresponds to the logarithmic bipartite fidelity (LBF) in Sec. 3. Interestingly the correction takes a simple 
form, that we test numerically with high precision in free fermionic systems. We also consider the effect of 
boundary changing operators, and provide an exact calculation of this term for an XX chain divided into 
two equal halves. Sec. 4 deals with the study of the entanglement entropy and the Loschmidt echo following 
a local quench. The corresponding imaginary time domain has the shape of "double pants". There is still 
a L _1 logL correction in this case, but intriguingly, it disappears after analytic continuation to real time. 
However, we show that a slight modification of the Loschmidt echo — dubbed "detector" in [30] — preserves 
the correction. We finish with some concluding remarks in Sec. 5. 

Some additional information is provided in three long appendices. Appendix A summarizes the 
connections between the bipartite fidelity, corner free energies, X-ray edge singularities, as well as some 
closely related recent attempts at measuring entanglement. The other two are more technical. In Appendix 
B we generalize the original calculation of the universal LBF finite-size function presented in [15], to account 

\ There may be more than one UV cutoff a in the different terms in the expansion, for simplicity here we look at the effect of 
a simple shift L — > L + a. It is also possible to perform the same analysis with more general power-law terms added to the free 
energy. 
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for possible changes in boundary conditions. Finally, we explain in Appendix C the exact calculation of the 
LBF for an XX chain cut into two halves, and how the asymptotic expansion allows to recover the CFT 
results. 

2. Criticality, corner free energy, and the extrapolation length in 2d statistical models 

In this section we focus on a generic statistical system defined on a lattice. We imagine that this system 
is at a critical point, such that the long-distance physics is captured by a Euclidean field theory that is 
conformally invariant in the bulk. The lattice spacing a will play the role of a UV cutoff. We are interested 
in a system with a boundary. 

2.1. The extrapolation length as a boundary perturbation by the stress-tensor 

In the study of critical phenomena, the concept of extrapolation length is used very often in the presence 
of a boundary in the system [25]. The idea is the following. Take a d-dimensional system on a lattice 
(aN) x (aZ) d_1 , with the plane (0, X2, £3, ■ ■ • , Xd) playing the role of the boundary. If the system is critical, 
it is described at large scale by some scale-invariant field theory in the bulk, with a scale-invariant boundary 
condition at the surface. In many cases, in order to relate lattice observables to their field theory counterpart, 
one needs to make the assumption that, in the continuum, the boundary of the system is at position x\ = — £ 
rather than at x\ = 0. The length £ > is dubbed "extrapolation length"; it is usually of the order of the 
lattice spacing a, and it plays the role of a UV cutoff when one deals with quantities that are sensitive to 
the presence of the surface. 




u 

For instance, consider the expectation value of some local observable on the lattice at position (xi, 0, . . . , 0). 
When x\ 3> a, its behavior is captured by a one-point function (<p(xi, 0, . . . , 0)) of the corresponding local 
operator (f> in the field theory. This one-point function is fixed by scale-invariance, and by the value of the 
extrapolation length £: (tfi(xi,0, . . . , 0)) cx l/\xj + £| , where h is the scaling dimension of the operator </>. 
Since this holds in the regime x\ ^> a ~ £, (<p(xi,0. . . ,0)) ~ l/|xi|' l (l — h£/x\ + ■■•), so the presence of 
the extrapolation length £ > determines the form of the leading correction to scaling for this one-point 
function. 

In general, lattice effects may be tackled by adding perturbations by local operators to the field theory 
action S. The perturbations may live in the bulk or at the surface of the system. One possible perturbing 
operator at the surface is the stress-tensor T^ v . Of course, other perturbations can appear at the surface, 
and should be taken into account, but let us start with the stress-tensor only. The perturbation by the 
component T±± = Tu along the d — 1-dimensional boundary can be written as 

S — > S + ^ J d^xn T ±1 _ . (4) 

The coupling £ has the dimension of a length; it is nothing but the extrapolation length itself. To see this, 
recall that the definition of the stress-tensor is that, under an infinitesimal transformation Xn 1— > x^ + £^{x), 
the variation of the action is 

SS = -i- J d d *T^s v , (5) 

where the factor 2n is a normalization convention (it is the standard convention in CFT). Now consider a 
transformation that moves the boundary of the system from x\ — down to x\ — — £, for instance S\(x) = £ 
for X\ < and e\ — for X\ > (and £2(x) = everywhere). Then the only non- vanishing component of the 
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tensor d^e v is d\£i(x) — — £S(xi), and it gives back the expression (4) when it is inserted in (5). Thus, the 
appearance of the extrapolation length £ is due to a perturbation by the stress-tensor along the boundary §. 

It is natural to ask what other boundary perturbations may occur. In general, this depends on the 
system, and on the local operators that are available in the field theory. First, one has to identify which 
conformal boundary condition describes the system in the scaling limit. Second, only operators that are 
irrelevant along the surface (in some cases marginal ones could also be present) can appear as perturbations, 
since relevant operators would drive the system towards a different conformal boundary condition under 
the Renormalization Group (RG) flow. Then the determination of the possible surface perturbations boils 
down to a standard RG analysis based on symmetry criteria and comparison of scaling dimensions of the 
operators present in the field theory. In this paper we will be mostly interested in the subleading corrections 
due to the stress-tensor; to study these corrections, we may assume that we are in a situation such that the 
stress-tensor is the least irrelevant boundary perturbation. 



2.2. Leading order of the corner free energy at criticality: the Cardy-Peschel formula 

From now on, we restrict our discussion to the two-dimensional case. In two dimensions, as pointed out by 
Cardy and Peschel [17], there is a universal contribution to the free energy of a critical system in a domain 
with a corner with internal angle 9. For completeness, in this section we quickly recall their result and its 
derivation; this is independent from the extrapolation length and boundary perturbations described in the 
previous section. The expert reader may skip this discussion and go directly to the next section, which 
contains new results about the extrapolation length and subleading corrections to the corner free energy. 




Figure 1. Sharp corner with internal angle 9. L is a typical length in the system. 



The Cardy-Peschel term can be computed as follows. The mapping z i— > w(z) 



maps the upper 



half-plane onto the corner shown in Fig 1. The holomorphic component of the stress-tensor in the corner is 

-2 



T(w) 



dw 
dz 
z 2(i-e/7r) 



T ( z ) ~ i^{w,z} 
T(z)- 
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(6»/vr) 2 f K ~' 12 2z 2 
Here we have used the transformation law of the stress tensor [31], 

-2 



T(w) 



dw\ 
dz J 



T ( z ) ~ ^{ w ' z } 



and the fact that (T(z)} = in the upper-half plane, {w, z} denotes the Schwarzian derivative: 



{w,z} = 



d 3 w/dz 3 3 / d 2 w/dz 



dw /dz 



dw/dz 



(6) 



(7) 



(8) 



The generator of dilatations in the corner is then given by the following integral over the dashed line in Fig. 1 
1 
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2ni 



wT(w)dw + c.c 



§ We are grateful to Nick Read for making this simple but crucial observation. This equivalence between the idea of an 
extrapolation length and the perturbation by the stress-tensor along the boundary plays an important role in [27]. 
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where C is a contour encircling the origin in the complex plane (we have used T(z) = T(z) on the real axis, 
and the standard boundary CFT trick that the non-chiral theory on the half-plane is like the chiral theory 
on the plane [32]). Now, the free energy F corncr scales as 



7 - I T-° 

-^corner \ J -' 



When there is no operator at the corner, ( Lq ) — 0, and one finds the Cardy-Peschel term 



F r , 



c 

24 Vtt 



logL. 



(10) 



(11) 



When the boundary condition is different on the two sides of the corner, one must insert a boundary condition 
changing operator [33] at the corner, which is a primary operator with a conformal dimension hb cc . In that 

case, (I/q) — h-b cc , and the Cardy-Peschel formula becomes 



t>bcc 



C 

24 V it 



logL. 



(12) 



2.3. Extrapolation length in the corner 

The question is now: what kind of subleading corrections to the corner free energy should we expect? Can 
we adapt the idea of an extrapolation length to a domain with a corner? In the picture shown in Fig. 2, this 
does not look completely obvious: the effect of the extrapolation length should be a shift of the sides of the 
domain, but it is not clear what one should do close to the corner. The question makes more sense when it 




Figure 2. How can one adapt the idea of the extrapolation length to a domain with a corner of angle 91 

is formulated in terms of the perturbation by the stress-tensor along the surface, 

S — > S + J T ±± dl, (13) 

which gives a correction to the corner free energy -Fcomcr - > ^corner + ^corner- Here dl is a line element along 
each of the sides of the domain. We want to evaluate the correction at the first order in £ for the corner 
shown in Fig. 2, 

^corner = ~~ / d W (T(w))-e t9 ^ f dw (T( W )) + C.C. (14) 



Here C\ is the horizontal side of the corner in Fig. 2 and C2 is the image of C\ under the rotation of angle 
9 (the other side on the figure). The factor e 10 comes from the Jacobian of that rotation. As we shall see 
shortly, one needs to be careful with UV and IR cutoffs in the integrals. But first, let us evaluate (T(w)). 
Again, it can be computed using a conformal mapping from the upper half-plane to the domain with the 
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corner, and using the transformation law of the stress-tensor. However, unlike the Cardy-Peschel term, one 
expects the subleading corrections to be sensitive to the full geometry of the domain, not only to the sum of 
the different angles. Therefore, we need to work with a slightly more precise form of the conformal mapping 
from the upper half-plane to the domain with the corner shown in Fig. 2. More precisely, we are interested 
in a generic simply connected dom ain with a corner of internal angle 9. We take a conformal mapping 
z i Y w(z) from the upper half-plane to this generic domain with a corner located at u>(0). Such a mapping 
always exists, thanks to the Riemann mapping theorem [34]. Close to the corner, the mapping takes the 
form: 

w(z)-w(0) = z e ^ f(z), (15) 

where f(z) is a function that is analytic around z = 0. Using the SL(2,M.) reparametrizations of the upper 
half-plane 1 1 , f(z) can be chosen such that it has the following expansion around z = 0: 

f(z) = l + f 2 z 2 + f 3 z 3 + ... (16) 

Such a function f(z) ^ est appears as soon as the domain is not exactly a single infinite corner 
{w € C; argz € [0,0]}, but differs from it at some scale L. For instance, one can think of some polygon 
with sides that have a length of order L. The typical scale L appears in the mapping z i-> z 8 ^ f(z) as the 
coefficient f2- typically, we have L ~ |/2|~ 27r/ ^- 

Now, with the mapping (15,16), we can compute the expectation value of the component of the stress- 
tensor (T(w)}. Again, we use the transformation law (7) to relate (T(w)) to the Schwarzian derivative 
of the conformal mapping z i— > w(z), and to the expectation value of the stress-tensor in the upper half- 
plane, (T(z)), which is zero (here we assume that there is no boundary condition changing operator on the 
boundary). After some simple algebra, we find that the correction to the free energy takes the form 

£c f~ Zuv f dw\~ 1 £c f Zia fdw\~ 1 

5F mme . = — dzl- {w ,z} + — dzl- {w,z} 



-zuv 

ZlR 



B | f C 

-,1+6/tt ~ J2 z 6/tt- 











PZIR 


1 + 




dz - 


r- £ 












'zuv 



z l+8/ir ~ J2 z 6>/tt- 



dz. (17) 



Here, each of the two integrals corresponds to one side of the corner. On each side, the variable w is integrated 
between an UV cutoff of order ~ a (lattice spacing) and an IR cutoff which is the system size, of order ~ L. 
The values of zjjv & n d zm are the inverse images of these two cutoffs under the mapping z t— > w(z); we 
are thus integrating over the variable z between zjjy > and zm > zjjv for one side of the corner, and 
between — Zjr and — z\jy for the other side. Note that, since w ~ z 6 ^, the UV and IR bounds for the 
variable z are respectively of order zuv ~ aV and zm ~ ~L~s . The precise values of Zm, zjjv depend on the 
details of the system. Here, what is important is that the first term in (17) gives a contribution of the form 
A + B'^/L where A, B' (and B and C) are numerical constants of order 0(1). The second term also gives 
such a contribution, unless 6 — This is the reason why = 2tt is so special: when 9 < 2ir, the corrections 
to the free energy are algebraic, with coefficients that are some mixture of the different cutoffs, and that 
are not particularly meaningful. As emphasized in the introduction, the case when there is a logarithmic 
correction is much more interesting; this happens when 9 = 2ir. 

2.4. Logarithmic correction in the case 9 — 2~k 

When 9 = 2ir, we see that the second term in (17) gives a logarithm piece when one integrates from an UV 
cutoff of order ~ = ^/a to an IR cutoff of order ~ = \[L. A proper evaluation of the Schwarzian 
derivative {w, z} of the mapping (16) leads to the following term: 



^-^Gorner A 



12tt./ .at 2 



= A + ^Jl\ og L/a + 0{l/L). (18) 

167T 

|| These are the special conformal transformations f(z) = °^t d , with a,b,c,d £ R and ad — be = 1. They leave the upper 
half-plane invariant, so if z 1— > w(z) maps the upper-half plane to the desired geometry, z 1— > tt>(£ (2)) also does. 
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The first term A is a non-universal constant of order 0(1) that depends on a and £. The second term is of 
order 0(L _1 logL); it is the one we are interested in. It depends on a single length scale, the extrapolation 
length £. In that sense, it is "semi- universal" : once £ is known, then this term is fixed entirely by the 
geometry of the sample, which is encoded in the single parameter / 2 (more precisely, L x / 2 , which is a 
dimensionless parameter). Note that such a statement does not hold, for instance, for the correction of 
order O(l), that is typically a ratio of several short length scales like £ and a. Eq. (18) is one of the central 
analytical results of the paper; we will show later that this logarithmic term is calculable in lattice systems, 
and we will find that it agrees remarkably well with numerics. 

To conclude this section, let us recast our result (18) in a more general form, valid for any simply 
connected domain with several corner of angle 2tt (we will later consider specific examples of such domains) . 
This more general form can be obtained by relaxing the condition (16), and by repeating the same arguments 
as above. The expressions involved in the calculation are slightly more complicated, but the arguments 
remain identical. For a given conformal mapping z i— > w(z) from the upper half-plane to the domain with 
the corners, we find the following correction to the free energy: 

SF COIaeI = A' + c f geom . |log(i/a) + 0(1/L). (19) 

The constant A' is of order O(l); the interesting term is the second one. The geometric factor is a 
dimensionless parameter defined by 



L 



24tt ^ 

2n corners 



(20) 



The sum runs over the different corners with internal angles 27r, which are located at the positions w(z c ), 
with the z c 's along the boundary of the upper half-plane. The phase factor e la " appears after some algebra 
when one relaxes the condition (16); a c is defined as 

) • (21) 




a c = arg 

We will discuss several examples of such geometric factors for specific domains below. 
2.5. Correction to the one-point junction 

The correction S — > S + 5S due to the extrapolation length (the stress-tensor) must also affect the one-point 
function of a primary operator (j)(w, w) close to the corner, 

(4>(w, W)) 5S = (4,(w)) + S (tfrw, w)) . (22) 



Figure 3. How does the extrapolation length in the corner modify the behaviour of the one-point function? 

The one-point function without perturbation (<p(w,w)) can be computed using the inverse of the conformal 
mapping z > w(z) [31]. In this section we compute the correction S (<f>(w, w)} at the first order in the 
extrapolation length £. The variation of the correlation function is, at the first order in JS 1 , 

S (<t>(w, w)) = ^ W,W) * ^° = - (SS <p(w, w)) + (SS) (<t>(w, W)) Q . (23) 
\ e /o 
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From now on, all the expectation values that will appear are taken with the unperturbed action S, so we 
drop the subscript '0' in ( . ) - The terms on the right-hand-side involve 5S, which is a sum of two integrals 
along the two sides of the corner, C\ and C2. Let us focus on the first side C\. Again, we have to integrate 
over the variable w from the UV cutoff (~ a) to the IR cutoff (~ L), and we use a conformal mapping from 
the upper half-plane to our domain, z \— > w(z), with w(z) ~ z®^ for z — > 0: 

-L 



8{4>{w,w)) Cl 



2m 
2tt 



dw' [{T(w')(j)(w,w)) - (T(w')) ((j)(w,w))] 



c.c. 



dw 

dz 



-2h 



dz' 



dw' 
~dz' 



(T{z')^z,z)) 



(24) 



In the second line, we have used the transformation law of the stress-tensor (7), the transformation law of 
the primary operator (f>(w,w) [31], 

-2I1 



4>(w, w) 



dw 



dz 



(z,z), 



(25) 



and the fact that the terms with the Schwarzian derivative cancel in the difference (Tip) — (T) (0). The 
contribution of the other side of the corner C2 is similar. Summing these two contributions, we see that the 
correction to the one-point function must be of the following form: 



8 (4>(w, w)} 
(<f>(w,w)) 



2tt 



dz z 



rji-e/n (T(z')4>(z,z)} 



2tt 



dz z 



>jl-G/it (T(z')(j)(z,z)) 



(4>(z,z)) 



c.c. 



(26) 



Like in the case of the free energy, if 6 < 2ir, this will ultimately lead to algebraic corrections of order 0(1), 
0(1/L), and so on. When 9 = 2tt, however, things are again more interesting, in the sense that we find a 
logarithmic divergence. So, from now on, we focus again on the case 8 = 2ir. We notice that (T(z')<j>(z,z)} 
is analytic at z' = (since z is not at the origin), and we expand 

O(z'), (27) 



dw' 

Hz 1 



1 

2z! 



so we find a logarithmic divergence coming from the pole at z' = 0. We have 
5 (4>(w,w)) 



(<j)(w,w)} 



2tt 



/L 1 (T(0)ct>(z,z)) 
dz 



= ^ log(£/«) x 



2z> (4>(z,z)) 
(T(0)<f>(z,z)) 



2n 



, 1 {T(0)<Kz,z)) 
,_^- L dZ 2z< (fan)) +aC - 



+ c.c + 0{l/L). 



(<P(z,z)) 

To evaluate (T(Q)cj)(z, z)) / (<j)(z,z)) we use the conformal Ward identity [31], 



(T(O)0(z,z)) - 



h 1 

"2 ~ 
• - 



z 



{<i>{z,z)) > 



(28) 



(29) 



together with the known form of the one-point function in the upper half-plane: (<j>(z, z)) oc l/(z — z) 2h . 
This leads to 



(T(O)0(z,z)) 



h 

72 



2z 



1 



Adding this expression and its complex conjugate, we find that the correction S (<f>(w , w)) due to the stress- 
tensor at the boundary is of the form: 

5((j>(w,w)} £/i(Imz) 2 



(<p(w,w)) 



log(L/a) + 0(l/L). 



(30) 



Note that, since 



L, the logarithmic term that we find here is of order \og(L/a)/L. Formula (30) is 



our second main analytic result. It shows that the extrapolation length generates a logarithmic correction 
to the one-point function in a domain with a corner of angle 2n. If the domain has several such angles, then 
the total correction is simply the sum of the ones associated to each corner. We will provide checks of the 
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formula (30) in Sec. 3.4. Let us conclude here with the following remark: our assumption in this section was 
that the conformal mapping z i— > w{z) form the upper half-plane has the expansion close to the corner with 
angle 2ir: w{z) ~ z 2 + . . .. This form is invariant under a composition with a special conformal mapping 
of the form z i-> z/{bz + 1), with b G M. This means that our result (30) must also be invariant under such 
a reparametrization. This does not seem obvious in (30); however one can check that Im z/\z\ 2 is indeed 
invariant under the reparametrization z t— > z/(bz + 1), b G K. 

3. Application to the bipartite fidelity and related quantities 

In this section we discuss the consequences of the general results of Sec. 2 on the bipartite fidelity of Id spin 
chains [15]. This quantity is an overlap defined with respect to a bipartition of the system. For a given 
bipartition of the sites of a spin chain A U B, the bipartite fidelity is the overlap between the ground state 
\A U B) of the Hamiltonian of the full spin chain, and the ground state \A ® B) of the Hamiltonian where 
all the interactions between the subsystems A and B have been switched off. We will focus on the logarithm 
of this overlap, dubbed "bipartite logarithmic fidelity", or LBF: 

T A . B = -log \{A (JB\A®B}\ 2 . (31) 

Below, the bipartition of the sites of the chain is chosen as follows: A corresponds to the first I sites, and B 
to the L — £ remaining sites (see Fig. 4). \A U B) is then the ground state of the total chain of length L, and 
| A® B) = \A) <g> \B), where \A) (resp. \B)) is the ground state of the chain A (resp. B) of length £ (resp. 
L — £). Some of our motivations for the study of the LBF are explained in Appendix A; it is a quantity 
which can be studied analytically by various methods, including CFT methods [15] and lattice integrability 
techniques [35, 36]. 

We start with the most generic case (Sec. 3.1), before generalizing to possible changes in boundary 
conditions (Sec. 3.2). The system considered here will be an open chain cut into two parts, as shown in 
Fig. 4. 



A B 
I 1 1 

< > i > 

i L-t 



Figure 4. Open chain of length L cut into two subsystems A and B of respective length I and L — I. 



3.1. Shape dependence in the generic case 

Let us first recall some basic facts about the LBF, following our previous study [15]. In imaginary time, the 
fidelity we are interested in can be expressed as a ratio of partition functions 

\(AUB\A®B)\ 2 = } Za - b J . (32) 

The corresponding geometries are shown in Fig. 5. In terms of free energies / = — log Z, the result reads 
•F a,b — 2fA,B — /a®b — Iaub- (33) 

Each of these free energies contribute to a bulk free energy (oc L 2 ) and a line free energy (oc L), but these two 
contributions are canceled by the linear contribution (33). In the end the leading (logarithmic) contribution 
comes from the Cardy-Peschel 2n corner in Ja,b (see Fig. 5). 
Let us introduce the aspect ratio 

x=|, (34) 
which we keep fixed in the following. We expect the LBF to scale as 

?(x, L) = °- log L + f{x) + g(x) 1 ^ + 0(1/L). (35) 
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I 


e 




2tt 


L-l 



Ia.b = - log Z A ,_ 





I 




L-l 



log Z A0B 




Iaub = — log Zaub 



Figure 5. Free energies f\ g, Ja®b an d Jaub involved in the calculation of the LBF. The leading 
logarithmic contribution comes from the corner with angle 9 = 2ir in /a,s- 



The leading logarithmic can be obtained as an application of the Cardy-Peschel formula (12): it is therefore 
a signature of the underlying critical behaviour. Note also that when approaching criticality, this term is 
replaced by (c/8)logA [36], where A is the correlation length, in the regime L A a. f(x) is universal 
finite-size function, and has been computed in [15]. It is given by 

2\, , , /. „ 2 



/(*) 



c 

24 



2x- 1 



log(l 



l-2x 



\ogx 



est. 



x j - ■ \ 1 - 

We are, however, more interested in the subleading term proportional to L~ 1 \ogL. 
transformation from the upper half-plane is [15] 



( ^ L 
w(z) = — 

7T 



x log(z + 1) + (1 — x) log 



zx 
l-x 



1 



(36) 

The conformal 
(37) 



and has expansion 

w(z) — iL{\ — x) 



Lx 



-z 2 + 0(z 3 ) 



(38) 

Since the coefficient of z 2 



- 1 



(39) 



2tt(1 - x) 

around z = 0. It is then straightforward to apply the Eq. (20) derived in Sec. 
in (38) is negative, the phase factor is e la — — 1. We get the simple result 

£c / 1 

24 \x{l - x) 

£ is the extrapolation length. Notice that the precise forms of Eqs. (35,36,39) rely on the technical assumption 
that the (conformal) boundary condition be the same everywhere. This assumption will be relaxed in the 
following subsection (3.2). 

Let us now check this formula in the XX chain at half-filling, a simple critical model with central charge 
c = 1. In general £ is a non- universal quantity, it depends on the specifics of the lattice model. Here it can 
be obtained from the standard bosonization argument that for a chain of L sites between x = 1 and x = L, 
the boundary conditions on the free bosonic field ip have to be (see e.g. [37]) <p(0) = <p(L + 1) = 0. We 
therefore have £ = 1. As is well known, this model can be solved exactly using a mapping onto free fermions 
[38], and the overlap we are interested in can be expressed as a determinant: 

(40) 



F(x,L) = -2 log det M, 



K l<i,j<L/2 

The matrix elements are given in Sec. (Appendix C), we refer to it for the derivation. Using standard 
numerical routines, it is then straightforward to evaluate T numerically for large system sizes. The numerical 
results are shown in Fig. 6. 

In practice we extract the L~ l \ogL term by fitting the linear combination 

F(x,L)-±]DgL-f(x) (41) 

to the form a + logi + aii^ 1 for system sizes L = 512, 1024, 2048, 4096, 8192 and fixed x. The data 

agrees perfectly with the CFT prediction (39), confirming our analysis. Interestingly, the shape is convex- up, 
as opposed to the concave-down form of the constant term f(x). Let us finish with two remarks regarding 
the extrapolation length £ 
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I 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

x = £/L 

Figure 6. Numerical extraction of the log(L)/L term g(x) in the XX chain at half-filling. The result is 
compared with the CFT prediction of Eq. (39) with c = 1 and £ = 1. 

• There exist a known exact mapping from an XX chain of length L at half-filling to two critical Ising 
chains of length L/2 [39, 40]. As a consequence [15], Fy^(x,L) = 2F\ s i ng {x, L/2). This result agrees 
with the well known value of the central charge c = 1/2 and an extrapolation length £ = 1/2. 

• It is possible to change the extrapolation length, while keeping the central charge constant. To do so, one 
can change the hoping amplitude at the two boundaries of the XX chain: H = — tj (^ +1 Cj + h.cj , 

with e.g. tj = 1 for j = 2, 3, . . . , L — 2 and t\ = tL-i = t / 1. We also checked numerically that our 
results are compatible with this. 

3.2. The effect of boundary changing operators 

It is also interesting to consider a slightly more general LBF, defined as 

FcW = -log| a<5 ( J 4u J B|| J 4>^ ®\B) lS | 2 , (42) 

where we have imposed different conformal boundary conditions ad, a/3, j5 on the states \AU B), \A), \B), 
as is shown in Fig. 7. 

a P 7 S 
I 1 1 

< > < > 

Figure 7. Bipartition of the chain with different boundary conditions imposed everywhere. Notice that the 
change between f} and 7 directly affects the leading Cardy-Peschel term. 

The scaling form of the LBF is changed to 

F= (| + Vr) \ogL + f aPl5 {x) +g a p~ / s{x) h ^ +0(1/L), (43) 

where the leading term is only affected by the change in boundary condition from (3 to 7, and is obtained 
from the Cardy-Peschel formula, Eq. (12). It is possible to generalize the result established in Sec. 2 to 
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take into account such changes in boundary condition. The L 1 log L contribution we are interested in is 
modified to 



%Pis{x) = £ x 



hsa ~ k + ( k ~ hafs )l + ( k ~ ^) r 1 



(44) 



where hij is the dimension of the operator changing the boundary condition from i to j. Notice that also 
the result of [15] for fapysiz) is modified in a non-trivial way. The exact expression for f a f)ySi &s well as the 
derivation of both results, is given in Appendix B. 

Such changes in conformal boundary condition can be achieved by explicitly changing the boundary 
condition on the lattice (e.g. returning a spin), but are sometimes only apparent in the effective field theory 
description. We will consider the two cases in Sec. 3.2.1 and Sec. 3.2.2, to check the main result of this 
section, Eq. (44). 



3.2.1. Returning spins in the Ising chain Let us consider the first excited state of an open Ising chain in 
transverse field [38] . We use the convention where the spins are measured in the basis of the eigenstates of 
the transverse field tr|. As is well known, it is twofold degenerate with the spins at both ends pointing in 
opposite directions, af = ±1 = — a x L . Excited states therefore naturally change the boundary conditions. 
For example, one can study the following modified LBF, with a — 8 =t and /? = 7 =4-, 

,2 



In this case we have h 



a/3 



\og\ n (AU B\\A) 
- h 7 s = 



hfi and h, 



n 

5a ' 



-it 



1 



x(l — x) 



0, so that we get 
24 L 



(45) 



(46) 



For the Ising universality class, we have c = 1/2, h^i = 1/2 [41], and we have previously found £ = 1/2. We 
present numerical checks of this formula in Fig. 8, using the same method as in the previous section. 




Figure 8. Numerical extraction of the logL/L term in an Ising chain in transverse field with boundary 
conditions H h, and comparison with the CFT prediction of Eq. (46) 

The results are again in excellent agreement with our prediction. Notice that the boundary condition 
changing operators have a huge impact on the result, since they change in this case the sign of g(x). 
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3.2.2. Changing the filling fraction in a Luttinger liquid Let us now turn our attention to a slightly more 
subtle case, where the change in boundary condition is not obvious in the underlying lattice model, but 
rather emergent in the continuum limit. The prototype Hamiltonian we have in mind is that of an XXZ 
chain in magnetic field 

H = J2 Ow+i + i + A ^f + ( 47 ) 

The magnetization sector M — j- J2j °~j m which the ground state lies is a non-trivial function of the 
magnetic field h. We are interested here in situations where M/0 and therefore h 7^ 0. As has been noted 
in Refs [42, 43, 44], non trivial magnetization sectors may induce changes in boundary condition. To see 
this, it is most convenient to consider the special case A = 0, that is equivalent to free-fermions hoping on 
the line. The magnetization is related to the fermion density by M = (p— 1/2), and the ground-state energy 
is given by 

E (L) = - £ cos C^L) (48) 
m=l V ' / 

up to an unimportant constant. The Fermi momentum is kp = pir and the Fermi speed Vp = sin hp. The 
low-energy physics can be accessed [3, 4] using the Euler-Maclaurin formula 

E (L) = aL + /3- 7 ^-[c- 24h p ] + 0(l/L 2 ). (49) 

c = 1 is the central charge. h p is the dimension of a boundary changing operator, given here by 

It is proportional to the so-called "phase shift" emphasized in the bosonization literature. Its presence means 
that although the boson height field obeys a Dirichlet boundary condition at both ends of the chain, there is 
still a height difference between the two ends, that needs to be taken into account. For general A the phase 
shift is to our knowledge not known analytically (see however Ref. [45] for closely related calculations). We 
refer to [44] for a more thorough description. 

For the ground-state \A U B), this means one can set ip(0) = 0, and therefore f{L) = A(f. Similarly for 
\A) (resp. \B)) we have ip(0) — and <p{£~) — A93 (resp. — and ip(L) ~ A92). We have then to 

take into account four boundary changing operators with dimension h a p — /i J g 7 = h y s = hg a = h p in the 
calculation of the LBF. Using Eq. (44), we find that the L™ 1 logL term is modified to 

^) = ^( 2 ^-V)(^(lb)- 1 )- (51) 

We performed checks of this formula in the A = 0, case. The results are presented in Fig. 9, for two filling 
fractions p—l/i and p — 1/8. 

The numerics show once again a perfect agreement with our formulas. We did not perform any numerical 
checks for general A, but we expect our formulas to be valid, with the more general 

R is the compactification radius^], obtained by solving numerically the integral equation derived from the 
Bethe-ansatz [45, 46, 44, 47]. S/n is the phase shift [44]. 

Let us now come back to the orthogonality exponent. As we have seen the leading behavior is given 
by J- ~ (c/8 + hp) logL, so that the LBF can be used to measure h p , without having to know a priori the 
Fermi speed vp. However the (single interval) entanglement entropy is not sensitive to this, and always scale 
as S ~ (c/6)log£. One might wonder why this is case. Such a behavior can be traced back to the fact 
that the replicated geometry used to calculate the EE doesn't really have any boundary, rather twists that 



% At half-filling it is known exactly [45] , R = \J2 — (2/n) arccos A 
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0.25 




Figure 9. Numerical extraction of the log L/L term in the XX chain away from half filling, and comparison 
with the CFT prediction of Eq. (51). Left: Filling fraction p = 1/i.Right: Filling fraction p = 1/8. In this 
case the dimension of the boundary operator is large enough to inverse the shape. 



allow to jump from one cylinder to the next. The LBF involves at least an explicit boundary (the cut), and 
the possible appearance of h p is intimately connected to the fact that for c = 1 compact theories, Dirichlet 
stands for an (infinite) parameter family of conformal boundary conditions [48, 49]: in bosonization language 
the forward scattering is exactly marginal. For c < 1 rational theories the number of conformal boundary 
conditions is finite, and a natural "free " boundary condition should not generate any change in boundary 
condition at the cut + . 



3.3. An exact result for the symmetric cut 

For a symmetric cut x — IJL — 1/2, it was shown [50] that the determinant (C.28) can be computed 
explicitly. Let us note N = pL the number of fermions in the chain, 9k = kn/(L/2 + 1) and 4>i = 1%/(L + 1). 
Then the fidelity is given by 

J-=-21og|Avi, (53) 

with 

nfii sin2 &k Uf=i sin f*^) rii<fc<i<jv/2 ( cos d k ™ cos 9 t ) 2 U^ITkn ( cos 4>k - cos 0,) 

D N = " " = ' /? N/ n „ = = ■ (54) 

[(L/2 + l/2)(L/2 + 1)] N ' 2 Yf k L\ Ul, (cos 9 k - cos & ) 

From this formula it is possible to perform a systematic asymptotic expansion, and prove the CFT result 
at x — 1/2. The details of the expansion, as well as the derivation of Eqs. (53,54) are given in Appendix C. 
The final result reads 

oo 

F{L)=Y J (lk\ogL + p k )L- k . (55) 

fe=0 

The equality is meant in the sense of asymptotic series. The 7/c are given by 

7o = ^ + h p , (56) 

7 fc -(-l) fe+1 (^-v)(2' £+1 -l) , fc>l. (57) 

7o is nothing but the Cardy-Peschel term, and 71 is the L _1 logL term. Both match perfectly the CFT 
prediction. The non universal [i k terms are much harder to obtain from the exact solution, and less interesting 

+ Except of course if one explicitly forces the matter, as we did in 3.2.1. 
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to us. We only computed fio in the appendix. It is quite striking that all the expressions for the 7& are 
very simple. Such an observation is consistent with the general idea that terms protected by a logarithm are 
less dependent on the numerous UV cutoff present in the systems. We did not attempt to calculate them 
using CFT for k > 2. The main reason for that is that they are slightly less robust than the L -1 logL; in 
particular they are not semi-universal, in the sense of Sec. 1.2. Note finally that this exact result also applies 
to the Ising chain, because Fxx{ x i L) = J r i S i ng (x, L/2), due to the exact mapping of Refs. [39, 40]. 

3.4- Correction to the one-point function: two examples 

In this section we provide numerical checks of our formula (30) for the 0(L _1 logL) correction to the one- 
point function close to a corner. The two geometries we consider are pictured in Fig. 10. 








•0 


(w, w) 






L 












L/2 

w = X 





(a) 







(b) 







Figure 10. There is a correction of order 0(L 1 logL) in the one-point function. Here are the two 
geometries we use as examples. 



3.4.1. One-point function in the plane with a slit Our first example is about the one-point function (4>(ui, w)} 
in the geometry shown in Fig. 10(a) . The domain is the complex plane minus a segment of length L: C\[0,L}. 
To compute the one-point function (<p(w,w)} in conformal field theory, we first need a mapping z v- > w(z) 
from the upper half-plane to C \ [0,L]. We chose the following one: 

w(z) = (58) 

L 

Using the inverse of this mapping, and the fact that the one-point function in the upper half-plane is 
proportional to l/|Im;z(u;)| 2 ' 1 , we find that the one-point function in C \ [0, L\ is, at the leading order: 

2/1 



(<f>(w,w)) 



dz 
dw 



2h 



(<t>{z,z)) oc 



w) 



(59) 



|f(l-f)| 5 ImyMiT 

Now we turn to the subleading correction of order 0(logL/L). The mapping (58) has the behaviour we 
want at small z, 

1 



w(z) ~ + ...), 



(60) 



so we can apply our formula (30) directly. We find that there is a subleading contribution 5 {4>{w, w)) coining 
from the corner at position w = 0: 



6{<t>(w,w)) _ £ /i Im \/f (1- f ) log(L/o) 



((j>(w,w)) 



\w/L\< 



0(1/L). 



(61) 
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There is another contribution coming from the corner at position w = L. This can be obtained simply by 
symmetry, by substituting w — > L — w. In total, the subleading correction to the one-point function coming 
from the two corners is: 

2 

1 1 



5 (4>{w, w)} £ h 

((j)(w,w)) 7T 




+ 



w/L\ 2 \l-w/L\ 2 



x l0g( ^ /a) + . . . (62) 



When w lies on the real axis, at a distance I > from the right corner, the results simplify into 

/ \ 2h 

(0) « -. =- .J L- 2 '\ (63) 



^ (^> , i \ iog(i/o) 



0(1/L). (64) 



It is possible to check this formula in the Ising chain in transverse field at the critical point, by studying the 
quantity 

(<TL+t) m = 2 (P L+ i) slit ~ 1, (65) 

where Pi = (1 + cr J :)/2 is the projector onto the state with up spin i. We measure the spins in the basis of 
the eigenstates of the transverse field Pauli matrices. Upon performing a Jordan- Wigner transformation, p 
projects onto to a state where the site i is occupied by a fermion. A lattice realization of the slit can be 
obtained by setting 

(n\p L+e nf=iPj\^) 

(P L +i} s nt = — , — 7 — — , 66) 

where is the ground-state of the chain. The sequence of L up spins is expected [29] to renormalize to a 
free conformal invariant boundary condition. Note that the denominator in (66) is nothing but the emptiness 
formation probability [51]; we refer to Ref. [29] for a detailed CFT treatment. Both terms in the ratio (66) 
can be obtained as simple determinants, using free fermions techniques [51]. In general we have 

(si\f\p Xi \n) = det (^- + 7^ — TTx), (67) 

\ ill x, i / 1 < i>J .< n ^ 2 2Tr(x t - Xj + 1/2) J ' y ' 

which allows for fast numerical computations on very large system sizes. For a fixed aspect ratio £/L, we 
observe that the one point function scales as 

K + ,) slit - 2 - = ^ (l+^A) 1 ^) +0(1/L% (68) 

so that the leading term is proportional to L^ 1 . By direct comparison with (63) we recover the dimension 
h = 1/2 of the fermionic field. The numerical results for f(£/L) are shown in Fig. 11 on the left, for L = 4096. 
We observe that the conformal scaling of (63) is perfectly obeyed. We also notice that the proportionality 
coefficient is very close to 1/4, a result that is probably exact. The identification of the subleading term 
requires some additional care. In practice, we fit (o'L+e)^ — § to the form a\L~ x + b 2 L~ 2 log L + a 2 L~ 2 for 
very large system sizes, and plot the ratio 62/^1- The data, shown in Fig. 11 on the right matches perfectly 
the CFT calculation (64) with the expected values h = 1/2 and £ = 1/2. 
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Figure 11. Left: Leading term in the one-point function. Right: Subleading L 1 log L term in the one-point 
function. In both cases the numerics match perfectly the CFT calculation. 



3.4-2. One-point function in the 'pants' geometry Our second example is the one corresponding to the 
domain shown in Fig. 10(b). A conformal mapping from the half-plane to this "pants geometry" is 



L ( 2ttz 2 
„,(*) = --log 1- — 



(69) 



With the inverse of this mapping, we get the leading order of the one-point function, exactly as we did 
previously for the slit: 

2h 

1 



((j)(w,w)) oc 



L |sinh=f| 5 



ImVT 



(70) 



Again, we have chosen the mapping such that it has the following expansion for small z: 

w(z) ~ z 2 (l + ^z 4 + ...), (71) 
so we can apply directly our formula (30) for the subleading correction due to the stress-tensor. We find 

S(<t>(w,w)) £h( 5: Rfi „, ImVT 



(<f>(w,w)) 




|sinh^| 



log(^/a) 
L 



0(1/L) 



(72) 



We now proceed to check these results numerically in the Ising chain for w on the imaginary axis, at a 
distance y > from the the corner (see Fig. 10). With h = 1/2 and £ = 1/2, the above formulas simplify 
into 



(<f)(w,w)) oc 
S(4>(w,w)) _ 



Lsin^cos(f + X 



1 



1 



sm 



EM. 



- 1 



L 



0(1/L). 



(73) 



(74) 



(4>(w,w)) u v ,.u L 

The results are shown in Fig. 12. We use a similar procedure to extract the contributions to the one-point 
function, but finite-size effects turn out to be much stronger, especially for the L^ 1 logL term. To remove 
them, we performed fits with more corrections, of the form a\L~ Y +b2L~ 2 log L+a2-L~ 2 +&3L~ 3 \ogL+asL~ 3 . 
The agreement with CFT (red) curves becomes then excellent. Notice that this particular form of correction 
is quite natural in such setups, as it is the same as for the exact calculation presented in Sec. 3.3. 
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Figure 12. Left: Leading term in the one-point function, and comparison with CFT (red curve). Right: 
Subleading L~ 1 logL term, and comparison with CFT (red curve). The convergence is better if one includes 
more subleading terms in the fit (black stars). 



4. Time evolution following a local quench 

This section is devoted to the study of time-evolution problems, following a local "cut and glue" quench. As 
before, we are looking for global observables that exhibit corner-singularities with angle 2ir. We will study 
here the Loschmidt echo and the entanglement entropy, looking for subleading corrections to the results we 
obtained in Ref. [30]. It will turn out that the application of our main result is more tricky here, due to 
interesting subtleties involving the analytic continuation to real time. 

Let us consider a quantum system prepared in the ground-state \A® B) of the chain cut into two parts. 
At time t = + the two subsystems are joined, so that the initial wave function evolves with the Hamiltonian 
Haub- 

\iP(t)) = e itHAUB \A®B). (75) 

It is now well established that such quenches can be interpreted using a quasiparticle picture [26]. For the 
local quench these are emitted from the region near the cut between A and B, and propagate at the Fermi 
speed vf [52, 53, 54]. They can typically be observed studying correlation functions, but we are interested 
here in global observables, such as wave function overlaps or the entanglement entropy. In a finite system, 
the quasiparticles bounce back on the boundaries, and all observables should become periodic in time (see 
Fig. 13). 



T[t/L) S(l/L) V{t/L) 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

t/L t/L t/L 



Figure 13. Time evolution of several observables following a local cut and glue quench. Left: Logarithmic 
Loschmidt echo F(t) for L = 128. Middle: entanglement entropy S(t) for L = 128. Right: Modified detector 
Loschmidt echo T>(t), slightly anticipating on Eq. (87), for L = 256. Each is periodic with period L/vp, as 
a consequence of the quasiparticle interpretation. We have set vp = 1 in the plots. 

Such behavior has been studied in detail in Ref. [30], for the entanglement entropy (EE) and the logarithmic 
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Loschmidt echo (LLE). Recall that the Renyi entanglement entropy (REE) is defined as 
1 



Sn = 



log tr A p" 



(76) 



1 — 71 

The von Neumann EE is the limit n — > 1, S\ = — tr^/dog/?. Note that the bipartition for the REE coincides 
with the position of the quench here, this will always be the case in the following. The LLE is defined by 

J-(t) = -logKV>(0)|^)>| 2 ; (77) 

it is a natural time-dependent counterpart to the LBF. Interestingly, the leading terms of both quantities 
can be computed within CFT. We have for two subsystems of the same size 



S(t) 

Ht) 



log 



log 



L . nvpt 

— sin 

7T L 
L . TTVpt 

— sin 

7T L 



(78) 
(79) 



These results have been checked in free fermion systems in Refs. [55, 30], see Fig. 13 (left and middle) for 
an illustration. Recently, Eq. (78) has also been confirmed numerically in Luttinger liquids away from free 
fermions [56] . Note that the similarity between the two formula can be deceiving, as their calculation turns 
out to be quite different. In particular, the relation J-(t) — (3/4)5(4) does not survive if La ^ Lb- In this 
section we check for the presence of a L~ x logL correction. We also set, for convenience, vf = 1. 



4-1. The Loschmidt echo 

The computation is performed in imaginary time. We need to evaluate 



J-(r) = -log\(A®B\e 



-tH 



B) 



(80) 

= -log |cst xZ(r) | 2 , (81) 

where Z(t) is the partition function of a "double-pants" geometry with spacing r between the slits, as shown 
in Fig. 14(a). We assume La = Lb = L/2 for simplicity here. The constant in (81) acts as a normalization, 
and is usually fixed when coming back to real time. Since the geometry has two 2-7T corners, we expect the 
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Figure 14. Multiple pants geometries representing the (logarithmic) Loschmidt echo (a) and a modified 
"detector" version (b). In each case the horizontal distance between the two-slits is r. 



LLE in imaginary time to scale as 



F(t) - - log 



L . 7TT 

— sinh — 

7T L 



cs t + 9(t) 1 ^ + 0(1/L). 



(82) 



The leading term has been derived in [30], and gives back (79) after analytic continuation. g(r) can be 
obtained from our main result in Sec. 2. A possible conformal transformation from the upper-half plane is 



W{Z ) = ^arccoth 

7T \ 7TZ 



L T7T 

—s- - coth— 
L 



with the two corners located at w = and w 
£c 1 



-r. Applying (20) to the conformal mapping, we get 



7(r) 



4 tanh(7rr/L) 



(83) 



(84) 
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for the prefactor of the L 1 log L term. Now we perform the analytic continuation r — > it + e: 

£c sinh(27re/L) — i sin(27rf/L) . , 

9 (it If =T X w„ k 7TV- (85) 

yy ' 4 cosh(27re/L) - cos(27rf/L) v ; 

e is a UV cutoff of the order of a lattice spacing. It is important to notice that (i) only the real part of this 
expression contributes to the LLE (see (77)) and (ii) the result can only be valid for times much larger than 
the UV cutoff e. So, although g(r) was of order one in imaginary time, its contribution becomes of order 
L _1 when coming back to real time. This means that we get, in the end, a correction 

£cire 1 log£ 

2 X 1 - cos(27rf/L) X £2 
to the LLE. Although the term (86) we just calculated should in principle be there, it is very difficult to 
identify it, even in a simple free fermion lattice model. This is mainly due to the other fast oscillating 
subleading corrections (see e.g. Fig. 13), which would need to be subtracted*. 

It is possible to circumvent this difficulty, and show that the logarithmic term we are investigating does 
appear in real time. To do so, we study a modified version of the LLE, called "Detector" in [30]. It is defined 
as 

V(t) = - log | (A <g) B\e itH \C ®D®E)\ 2 , (87) 

where | C <g> D eg) E) is now the tensor product of three ground-states. We take the respective sizes of C, D 
and E to be L/4, L/2 and L/4. See Fig. 14(b) for the corresponding space-time geometry. This quantity 
has 3 corners with angle 27r, and we expect it to scale as 

V(t) = ^ log L + d(t)+ est + e(t) x ^|^ + 0(1/L). (88) 
8 L 

The constant term is given in the plateau region by [30] 

d(t) = ~ log(2cos(27rt/L)) + 21og (l + 2 cos(27rf/L) - 2\/2cos(7rf/L)\/cos(27rf/ J L)) ,(89) 

and agrees with the numerics (see Fig. 13, right). To compute the subleading correction we need the 
conformal transformation from the upper-half plane to the detector geometry. It is given by 

W (z) = ^[log(l-z 2 )-21og(l-c 2 z 2 )] , r =i |log(4 W V 2 -l)). (90) 

Adding the 3 contributions coming from the three corners in Eq. (20), we find 

e 2irr/L + 2 % /2e- 7rT / L ,/cosh 2 ^ 




" r >"l 2+ cosh^ I 

Performing the analytic continuation we get 

which is of order one, as was hoped. Eq. (92) is much easier to test numerically, because it is less subleading 
than the analog term for the Loschmidt echo, and the results are not obscured by additional fast oscillation, 
at least in the first "plateau" region < t < L/4 (see Fig. 13). In practice, we extract e(t) by studying the 
linear combination 

(V(t)-V(0)-[d(t)-d(0)})xJ-r. (93) 

logL 

The results are shown in Fig. 15 for an XX chain at half-filling. Eq. (93) should go to (92), but the time 
dependent L _1 contribution to the detector translates into a correction of order 1/logL in (93). This makes 
the convergence very slow, as observed in the figure. However, studying this quantity for various large system 
sizes at fixed t/L and once again performing fits to b-\ + /_i/logL allows to get much better results. These 
are the black crosses in the figure, and they agree very well with our CFT prediction. 

* Note that (86) is not even the leading time-dependent contribution coming from the stress-tensor. A more thorough calculation 
taking into account all the other non- logarithmic terms yields -F(t) = | log | & sin ^ | -\-cst + + £ e7r cot^t/L) 2 j Q g 

Such a contribution is clearly not the only one of order 1/L, as it has period L, and cannot describe the fast 0(1/ L) oscillations 
with period ~ 1 observed in Fig. 13 or [30]. 
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Figure 15. Extraction of the time-dependent logL/L term of the detector in the plateau region. Big 
colored dots are the linear combination of Eq. (93), for system sizes L = 128, 512, 2048, 8192. This converges 
very slowly to the CFT prediction (red thick line). However, performing a fit to b_i + f-i/ logL improves 
the results a lot, and yields a very good agreement with CFT. The results are shown in black stars, for a 
few values of the reduced time t/L. 



4-- 2. Entanglement entropy 

The subleading corrections to the entanglement entropy in the ground-state are now relatively well 
understood [57, 58, 59], and it is natural to ask what kind of corrections are present after the "cut and 
glue" quench. Here we look for possible corrections of order O (logL/L) or 0(logL/L 2 ), coming from the 
stress-tensor at the boundary. To compute these subleading contributions, we proceed as follows. We use 
the trick of Calabrese-Cardy [54], namely we observe that the one-point-function of the twist operator on 
the world-sheet behaves as the one-point function of a usual primary operator in CFT. More precisely, we 
have 

trp" oc ($„O,u0) , 
where <&„ is a primary operator with conformal dimension 



h„ = — 
12 



1 

n 

n 



We evaluate 



(<$> n (w,w)) = 



(A(g>B\e- HT2 $ n (w,w) e- HTl \ A(g> B) 



(94) 
(95) 

(96) 

(97) 
(98) 

Eq. (96) is nothing but the one point function in the double pants geometry shown in Fig. 14(a), with 
r = T\ + T2- The corresponding mapping z H> w(z) is given by Eq. (83), with the two corners at w — and 
w = —t. The inverse is 



(A <g> B\ e-^Cn+Tj) \ A ® B) 
in imaginary time, and only later will we perform the analytic continuation 
Ti — > e + it, 
t-2 — > e — it. 



z{w) 




(99) 



The leading order of the one-point function was computed in [30], and we wish to use our result in Sec. 2.5 for 
the subleading logarithmic correction. There is however an important subtlety that needs to be taken into 
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account: the one-point function of the operator <£>„ is calculated in a regime where the spacing r = t\ + t 2 
between the two slits will become very small after the analytic continuation. This is very different from the 
situation we considered before, and it is not clear whether vL should still be the natural IR cutoff: naively 
one would now expect s/r to play this role. 

This intuition can be confirmed by an explicit computation of all the contributions from the stress-tensor 
to the one-point function in this particular geometry. To do so we perturb on the intervals [e, A], where e 
is a cutoff of the order of the lattice spacing, and A 3> L (the other perturbation on [—A, — r — e] can be 
deduced by symmetry). The corresponding contribution to the one-point function is (see Sec. 2.5) 

z(e) (resp. z(A)) is the image of e (resp. A) through the inverse conformal mapping (99). Using the Ward 
identity, we get the integral of a rational fraction, which can be performed. Adding the contribution coming 
from the other slit we get, after long algebra: 
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log 
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(101) 



El 
L 



The first term is of order L~ l logL provided t/L is of order one. However, after the analytic continuation 
n — > s + it, T2 — > e — it, with e another cutoff of the order of a lattice spacing, we observe that it becomes 
of order L~ 2 . Taking into account all contributions, we get 

S ($n) . , c , (2t/L - 1) cot 

— — =CSt + 27T^„ — 

\™n) Li 

In the end, the correction to the entanglement entropy is, using S = — lim T , 

est' £cir(l-2t/L)cot 
~L + 



-O(L-) 



Sit) = - log 



L TTt 

— sin — 

7T L 



est 



L 



,i d n tip n , 
+ 0(L- 2 ). 



(102) 



(103) 



As for the Loschmidt echo, we get that the leading time-dependent correction coming from the stress-tensor 
is proportional to L . The only qualitative difference between the two quantities is that the logarithmic 
term has transformed even more dramatically under the analytic continuation, from L~~ l \ogL to L~ 2 . Once 
again, (103) is very difficult to test numerically, as it is not the only contribution of order 1/L. 

Of course, it is a priori not obvious that the analytic continuation method should give correct results, 
especially when studying subleading contributions. Although it induces several intriguing subtleties in the 
calculation of the Loschmidt echo as well as the entanglement entropy, we have found no evidence that it 
could give incorrect results in situations were the leading terms are universal. It should however be stressed 
that this approach, even with more perturbations by local operators, always gives subleading corrections 
with period L; it cannot explain the fast subleading oscillations [53, 30, 56] with period of order one, that 
can also be observed in Fig. 13. 



5. Conclusion 



In this paper, we have studied possible corrections to the corner free energy of two-dimensional classical 
systems at criticality. We have focused on the very generic perturbation by the stress tensor at the boundary. 
In general this perturbation produces corrections of order L~ l to the leading universal Cardy-Peschel term. 
However, when the angle is — 2ir, the leading contribution from the stress-tensor becomes of order log L. 
The prefactor is proportional to the (non universal) extrapolation length £, but also the central charge as 
well as a universal function solely determined by the geometry. We have argued that such a correction is 
much more interesting than the usual power- laws of the form L~ a : in a sense the logarithm decouples this 
term from the other subleading corrections. This L^ 1 logL term also appears as a correction to the one-point 
function of an operator near the corner. 
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We have then applied this result to overlaps of one-dimensional critical wave functions, such as the 
bipartite fidelity. This was our initial motivation for the study of such corner free energies. Interestingly, 
it is possible to compare our CFT calculations with numerical evaluations of the fidelity in free fermionic 
systems. We have found a remarkable agreement between the two. For the special case of a system cut into 
two subsystems of the same size, it is even possible to obtain the fidelity in closed form for the XX and Ising 
chains, using a variant of the Cauchy double alternant formula. Performing an asymptotic expansion allowed 
us to recover the CFT result in this case. We also tested our analytical formula for the one-point function 
of a primary operator near the corner, and found again a perfect agreement with lattice calculations. 

Finally, we also explored the consequences of our result on the time evolution of several global 
observables, following a local cut and glue quench. It turns out this case requires some care, because of 
the final analytic continuation to real time. We computed the logarithmic correction to the Loschmidt echo 
in imaginary time, but found that it becomes of order L~ 2 \ogL after analytic continuation. However, certain 
"detector" deformations of the Loschmidt echo preserve the L _1 logL term, and we managed to establish 
their presence in a lattice model, confirming the validity of our approach. For the entanglement entropy 
there is a priori an analogous contribution to the one point function, but we showed that it becomes of order 
L~ 2 in real time. 

There are several interesting directions for future research. First, it would be interesting to study the 
sub-leading terms of the form L~ k log L, k > 2 that we identified in the exact lattice calculations. Such 
terms appear natural if one thinks of the second order effect of T(z) on the free energy, but there should be 
some other contributions involved, as they are not semi-universal. A precise determination of the possible 
structure of logarithmic corrections within CFT is left as an interesting question. It is also possible to study 
other Id quantum observables that exhibit 2tt corner singularities in imaginary time. One such example is 
the emptiness formation probability [51] in critical spin chains described by minimal models, where we expect 
the presence of a similar logarithmic corrections. A detailed study of this will be presented elsewhere [29]. 
Another interesting problem would be to study the general structure of the subleading corrections following 
the cut-and-glue quench, to see whether or not the analytic continuation still holds, and determine which 
operators generate the leading terms. Finally, it would be desirable to numerically check our predictions in 
interacting systems. The other corrections would presumably take a more complicated form, but we predict 
that the L _1 logL should still be there. For example in an attractive Luttinger liquid, it should still be the 
leading correction to the free energy, and, as a consequence, to the bipartite fidelity. 
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Appendix A. Measures of bipartite entanglement, orthogonality catastrophes, corner free 
energies, and all that 

In this appendix, we gather some of our notes on the relations between corner free energies, the bipartite 
fidelity, X-ray edge singularities and the closely related topic of the measurement of bipartite entanglement; 
these were our initial motivations for the work presented in this paper. Strictly speaking, this appendix 
does not contain new results, and most of its content has appeared in papers by other authors, in particular 
Cardy in [60]. It may, however, still be useful to some readers. 

Recently, there has been some interest in measuring quantities that characterize the bipartite entanglement 
of many-body systems [60, 61, 62], such as the Renyi entanglement entropies. Recall that these entanglement 
entropies are defined with respect to a given bipartition of the full Hilbert space Ha ®Wb of a quantum 
system. In this note, the quantum system will always be a 1 + Id system, such as a spin chain or a 
tight-binding model of fermions on a line. The subsystem A will be a segment in this line, and B is the 
complementary subsystem. This spatial partition AL) B induces a bipartition of the Hilbert space Ha ®T-Lb- 
If the full system is in some given pure state \tp), then the subsystem A is described by the reduced density 
matrix pa — tr# \ip) (tpl, where the partial trace is taken over all the degrees of freedom that lie in the 
complementary subsystem B. The Renyi entropies are defined as 

S n = log(trpl) . (A.l) 

1 — n 

In what follows, n > 1 is an integer, although strictly speaking the Renyi entropies could also be defined for 
non-integer n; in particular, the limit n — > 1 gives the von Neumann entropy. 

The difficulty that arises when one tries to relate the Renyi entropies to physical observables comes 
from the n th power of the reduced density matrix. Physical observables should be linear in pa, of the form 
(O) = tr [p A 0], and not of the form tr [p n A 0\. 

Appendix A.l. Replica trick, and the swap operator 

A key observation, which is at the heart of most of the work in this field, is that tr p\ can be obtained as a 
linear expression involving the density matrix of a replicated system. To see this, imagine that we have n = 2 
replicas of the full system ("Hai 'SiW-bi) ® ('Ma2 ®T~Lb2), in the pure state \ip) <8> IV*)- The density matrix that 
we are going to use is p® 2 = \ip) (ip\ (g> \ip) (ip\. One introduces the "swap" operator (also called "switch" or 
"twist" or "permutation operator" in various references), noted Sa, which maps canonically Hai onto 1-La2 
and vice- versa, and acts as the identity oti'Hbi®T~Lb2- More explicitly, this means that 5,4 is linear, and it 
acts on the basis states as 

S A {\n) Al ®\m) m )® {\p) A2 ®\l) B2 ) = {\p) A1 ® M B1 )® (|">A2® lOua) • ( A - 2 ) 
The expectation value in the state \ip) <X> of the swap operator, (Sa) = tr [Sa • p® 2 ] , is equal to 

(S A ) = tr nA P% (A.3) 

where the right hand side stands for the initial (non-replicated) system. This observation is easily extended 
to any integer n > 1, for a replicated system with n copies. The swap operator, in this case, maps T-Lau onto 
T-LAk+i, where k labels the n copies (say from 1 to n, modulo n). It is thus possible to get all the Renyi 
entropies S n by computing/measuring the expectation values of these swap operators. This is a trick that 
is useful analytically (see, e.g., [10, 12, 63]) as well as numerically [64]. A natural question is whether the 
swap operator can be implemented physically in a realistic system. 

Appendix A. 2. Action of the swap operator on a local Hamiltonian 

At this point, however, it is not clear that the swap operator may find a reasonable implementation in a 
concrete physical system, because as it is presented above, it seems to be a highly non-local operator. The 
next step is to re-express (Sa) as an overlap between the ground states of two Hamiltonians that differ only 
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by a few local terms [60] . Let us consider a concrete model to illustrate the idea. Take a homogeneous open 
chain with L sites, populated by free fermions with nearest-neighbor hopping 



L-l 
i = l 



[cjci+i + ■ ( A - 4 ) 



We chose the subsystem A as the first L/2 sites (i — 1, . . . , L/2, and we assume that L is even). The 
Hamiltonian H is equal to Ha + Hb + H AB where Ha (Hb) is the Hamiltonian of the left (right) subsystem, 
and the interaction term between A and B is H AB — t(c L ; 2 c i/2+i + c l/2+i c £/ 2 )' ^ ne ^ ac ^ that H AB is 
local, namely that it couples degrees of freedom that are located close to the cut between A and B only, is 
important: as we shall see shortly, it allows to reformulate the expectation value of the swap operator as an 
overlap between the ground states of two Hamiltonians which differ only by local terms. First, we introduce 
the n replicas of the system; each replica comes with the Hamiltonian 



L 

H k =t^ 

i=l 



\c\ k c i+1 , h + 4 +1>k a, k (A.5) 



— H A k + H B k + H AkuB k k = l,...,n, 

and the Hamiltonian of the full system, acting on the n replicas, is simply 

n n 

flrepl. - J2 Hk = E [ H *° + H B k + H^k] ■ (A.6) 
fe=l fe=l 

The ground state of this Hamiltonian is IV'rcpi.) = ® • • • ® 1^) ■ Now let us conjugate iT re pl. by the swap 
operator. This defines the "swapped" Hamiltonian, 

^repl. = S A ■ #repl. ' S A , (A. 7) 

which has the ground state |^epi ) = ^ A IV'repl.)- The ground state expectation value of the swap operator 
may thus be viewed as the overlap 

(S A ) = (CpiMcpL). (A.8) 

Such an overlap is, by definition, the (square root of the) probability of observing the system in the ground 
state of H' repl if it is initially prepared in the ground state of i? rep i.. In the recent literature, such out-of- 
equilibrium problems, where one suddenly changes the Hamiltonian H — > H' , have been dubbed "quantum 
quenches" [26]. In the present case, the two Hamiltonians -ff rep i. and H' rc ^ differ only by a few local terms, 

n n 

Kcpl. = E ^ HAk + HBk ^ + f E [ C L/2,k+l C L/2+l,k + C L/2+l k C L/2 , k+1 , (A.9) 
k=l k=l 

so this is one example of a local "quantum quench" . It is convenient to picture the two systems described by 
these Hamiltonians as in Fig. Al. It is clear from the drawing that the difference between i? rep i. and H' Tepl 
is local in the sense that it involves only terms that are located in the "center" of the replicated system (see 
Fig. Al). Therefore, the problem of measuring the Renyi entropy S n has boiled down to measuring the 

overlap (ip' xep i VW-) in a "local quench" H Iep i -t H[ cpl . 

Appendix A. 3. A simpler version of the "swap" quench: the "cut- and- glue" quench 
Implementing the "swap" quench H lcp \ — > -ff^ op i requires a number of skills: 

(i) one needs to create n identical copies of the quantum system, 

(ii) one also wants to be able to switch on/off one given coupling at a specific location in the middle of the 
system, namely to cut or glue two given subsystems instantaneously 

(iii) it is not sufficient to be able to cut or glue two of the subsystems; one actually needs to cut/glue 
simultaneously n + n subsystems, coming from the n replicas. 
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1 2 




Figure Al. An illustration (for n = 3 replicas) of the "swap" quench discussed by Cardy [60]: the 
Hamiltonian H rcp i (left) is instantaneously changed to ii^cpl ( r 'S n t)- Each shaded area corresponds to one 
replica of the original spin chain/tight-binding model. Each black line corresponds to an interaction between 
two sites. The key observation is that the difference H' Icpl — H rcp i, involves only a few local terms in the 
center of the picture. 

We note that it has been claimed recently that this type of quantum quench may be realized in cold atoms 
systems [61, 62]. At this point, it is worth emphasizing that, despite the modern vocabulary attached to these 
problems in the recent literature, local quantum quenches actually have a long history; their study dates 
back to the 1960's, with such phenomena as the X-ray edge singularities and the Anderson orthogonality 
catastrophe. This motivates us to look at a simplified version of the above "swap" quench, which does not 
involve replicas, and whose connection with the older class of problems is more transparent. 

To get a "baby- version" of the local quench H' rcpl —> H xep i , we simply drop the requirements (i) and (iii) . 
We do not introduce replicas of the original system, but we require the ability (ii) to switch on/off suddenly 
all the interactions B.\ B between the two subsystems A and B. Then the quench without the replicas is 
simply 

H = H A + H B -> H' = H A +H B + H A B . (A.10) 

Again, this is a local quench, in the sense that H' — H = H A B involves only a few terms that are located 
close to the cut between A and B. Like in the more complicated case of the "swap" quench, one can be 
interested in measuring the overlap between the two ground states of H and H' . Let us refer to this local 
quench as the "cut-and-glue" quench; it is a problem that has been studied in recent papers [54, 53, 15, 30]. 
Let us emphasize that the overlap between the two ground states H and H' is not a Renyi entropy like in the 
case of the "swap" quench. It is, by definition, an overlap, dubbed "bipartite fidelity" in [15]. Although it is 
not a Renyi entropy, it is a quantity that is interesting in its own right, and it can be used to characterize 
the entanglement between the two subsystems A and B. Measuring this quantity should be simpler than 
measuring the Renyi entropies, since the "cut-and-glue" quench is simpler than the "swap" quench. 

Appendix A. 4- Orthogonality catastrophes 

The overlaps between the ground states of two Hamiltonians which differ by a local term is a topic that 
was pioneered by Anderson in a seminal paper [65]. He showed that the two ground states of a Fermi gas 
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Figure A2. The "cut-and-glue" quench may be viewed as "baby-version" of the quench presented in Fig. 
Al. It does not involve replicas, but only the ability to switch on/off instantaneously the interactions 
between A and B. The goal is to measure the overlap between the ground states of the Hamiltonians 
H ' = H A + H B + H AlJB (top) and H = H A + H B (bottom). 

in the presence/absence a localized impurity potential have an overlap which decays algebraically with the 
system size L. This is the "orthogonality catastrophe" in its original formulation for a fermionic system. 
The exponent of the algebraic decay x, 

IW>I <* (~)*, (A.11) 

depends only on the low-energy properties of the system. Here a is an UV cutoff (typically, the inverse 
bandwidth). For free fermion systems, x is given by the phase shift S caused by the impurity potential, 
evaluated at the Fermi momentum: x = | (S(kp) /^)" ■ The orthogonality catastrophe in systems with a 
Fermi surface has been studied extensively since the 1960s. It is not a phenomenon that is restricted to 
systems with Fermi statistics though; the algebraic decay of the overlap between \^f) and l^') may also be 
interpreted as a signature of the underlying universal critical behavior. After reduction to s-wave scattering 
of the fermions on the impurity potential, and keeping only the radial component of the wave packets, the 
system can be reformulated as a 1 + Id semi-infinite gapless system, with a boundary condition related to 
the impurity potential through the phase shift S(kp). The boundary condition is therefore different when 
the impurity potential is on/off. A precise analysis of these boundary conditions leads to the orthogonality 
catastrophe and to the exponent x. This approach has been developed with great success for a wide variety 
of quantum impurity problems. 

The "cut-and-glue" quench in a 1 + Id system belongs to this class of impurity problems; it corresponds 
to the ideal case when the impurity completely decouples the two subsytems. The "swap" quench may also 
be viewed as an impurity problem, but the effect of the impurity is more convoluted, as it switches the 
couplings between n replicas of the system. The important point here is that both the "cut-and-glue" and 
the "swap" quenches in gapless systems are associated with orthogonality catastrophe exponents, x cut / s \ UQ 
and a; swap .In either case, the overlap between the ground states before and after the quench H — > H' obey 
a similar power-law decay with the corresponding exponent x. In systems with an energy gap, there is a 
finite correlation length £, and the two ground states can differ significantly only in a region of size ~ £ 
localized around the the cut between the parts A and B. In that case, one expects |(\& \*f?') | oc (a/t;) x as 
one approaches the critical point. 

Appendix A. 5. Measurable quantities related to the overlaps 

Strictly speaking, the (squared) overlap \W) | 2 is a measurable quantity itself, as it is, by definition, the 
probability of observing a system — initially prepared in the ground state of H — in the ground state of the 
new Hamiltonian H' immediately after the quench H — > H' , Nevertheless, this probability is usually not 
directly accessible in realistic many-body systems, especially when they are gapless. In the latter systems, 
the orthogonality catastrophe exponent is usually more easily accessible than the overlap |(^ \ ^') | itself. 

Appendix A. 5.1. X-ray edge singularity, and absorption spectrum of an auxiliary two-level system: it is 
very well-known that the orthogonality catastrophe exponent appears as a power-law singularity in the X- 
ray absorption spectra of metals. In this problem, the presence/absence of a localized core electron can be 
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viewed as a two-level system. In the presence of the core electron, the electrons close to the Fermi surface 
are in the many-body ground state |?/>). An X-ray photon can remove the core electron and therefore create 
a localized hole. The Hamiltonian for the electrons close to the Fermi surface is then modified to incorporate 
the localized potential created by the hole. This is precisely a local quantum quench problem. Since we know 
that the corresponding orthogonality catastrophe exponent appears in the absorption spectrum of the to- 
level system, it is a very natural idea to generalize this setting to any sort of local quench H — > H' in gapless 
systems, including the "cut-and-glue" and the "swap" quenches. One can imagine that the evolution of the 
system is governed by a Hamiltonian H or H' , depending on the presence/absence of the core electron, which 
may be represented by a two-level system, with two states |0) and |1), and energies and Eq respectively. 
The total Hamiltonian for the full system ([gapless system] (g) [two- level system]) is 

H <E> |0) (0| + H' ® |1) (1| + E Q l® |1) (1| , (A.12) 

such that the gapless system evolves with a Hamiltonian H if the two-level system is in the state |0), and 
with the Hamiltonian H' if it is in the state |1). If the state of the two-level system changes suddenly, one 
gets a quantum quench H — >• H' . To be more concrete, we consider again the toy-model (A. 4), in the case 
of the "cut-and-glue" quench. We imagine that the hopping term between the sites i = L/2 and i = L/2+ 1 
is coupled to the two-level system, 

.fftot. = t ^2 [4^+1+4+1^ 

When the two- level system is in the state |0), the left and right subsystems are decoupled, while the chain 
is homogeneous when it is in the state |1). The correlation function of the two- level system is related to 
the overlap between the two ground states in the presence/ absence of the impurity potential V, or, more 
precisely, to the orthogonality catastrophe exponent. For critical systems with a dynamic exponent z = 1, 
on time scales such that a <C vt <C L (v is the velocity of the massless excitations, a is an UV cutoff and L 
is the typical size of the system), one can show that 

(S-(t)S + (0)) = (-01 e- iHtat t \f) oc e- lEot x t~ 2x , (A.14) 

where S + |0) = |1), S~ |1) = |0) and S + |1) = S~ |0) = 0. The important point here is that the exponent x 
in (A.14) is the same exponent as in the orthogonality catastrophe above. The absorption spectrum of the 
two-level system thus has a power-law singularity from which the exponent x can be extracted: 

J dt e* Et (S-(t)S + (0)) ex \ E _l oV _ 2x - (A-15) 

Appendix A. 5. 2. Rabi oscillations of the two-level system: The following setting was discussed recently in 
[61]; it is essentially the same as the one in the previous paragraph, but in a different regime. Consider first 
the two-level system alone. Imagine that one switches on a small tunneling term T S + + T* S~ between the 
two states of the two-level system. If the system is initially (t = 0) prepared in the state |0), then at time 
t > 0, the probability of observing it in the state |1) is 

p\o)^\d (*) = {Eo/ ly 2 + w sin2 [* "/ 2 ] « = VWF+WW 2 (A.i6) 

Now, consider again the total Hamiltonian (A. 13), and assume that the tunneling amplitude |T| and the 
shift of the ground state energy Eq are both much smaller than the energy gaps of H and H' . In that limit 
(and only in that limit), one can project the system onto the two-dimensional subspace generated by |0) (g> 
and |1) \ip') . The tunneling term, when it is projected onto this two-dimensional subspace, becomes: 

V suhsp . -[TS++ T* S-] ■ V suhsp . = T W |V) S+ + (T (tf )* S- . (A.17) 

Thus, in the presence of the chain, the population of the two- level system oscillates as in (A. 16), but with the 
modification Eq — > Eq and T — > T (ij)' \ tp) . In principle, this can be used to extract the value of the overlap 
\(ip' \ip) |; if this can be achieved, it is more interesting than being able to measure the exponent x only, 
since the overlap itself is a more refined quantity than the exponent. However, notice that the restriction of 
a many-body system to a two-level system is quite a crude approximation, especially in a gapless system, 
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where the energy gap goes as 0(v/L). In particular, the shift of the ground state energy between H and H' 
may be hard to control with such accuracy (for instance, in the case of the "cut-and-glue" quench, the main 
contribution to the shift of the ground state energy comes from the change of boundary conditions at th c 
ends of the subsystems A and B, and this is of higher order; there is also a change of the Casimir energy 
which is itself of order 0(v/L); for the "swap" quench, although such contributions would ideally cancel, 
they may remain if the n copies of the system not perfectly identical) . 

Appendix A. 6. RG discussion of the perturbation H — > H' 

One-dimensional gapless systems with a dynamic exponent z = 1 are well described at low energies by field 
theories that are massless and Lorentz invariant (with the velocity of the low-energy excitations, v, which 
plays the role of the speed of light). It is most convenient to think of the original Id system as an equal-time 
slice of a 1 + 1-dimensional world-sheet. Both the "cut-and-glue" quench and the "swap" quench possess a 
nice interpretation in terms of the geometry of the world-sheet. 



In the field theory language, the difference between the Hamiltonians H' and H corresponds to a local 
operator, which we note A <f>(x) (<fi(x) is some operator in the field theory, and A is a coefficient). This 
operator perturbs the world-sheet action along a line corresponding to a fixed position x = x . In cuclidcan 
space-time, the action of the theory becomes 



This means that, while the original action S describes a homogeneous system on the euclidean world- 
sheet, there is now a defect along the line x = x - One needs to discuss whether the operator 4>{x,t) is 
relevant /irrelevant /marginal in the RG sense, when it is integrated along a line. If 4>(x,t) is irrelevant along 
the line, then it cannot induce a significant change of the ground state at large distance, therefore there is 
no orthogonality catastrophe. On the contrary, when (/>(x,t) is relevant, then the system flows towards a 
new IR fixed point, and the exponent x must be calculated at this IR fixed point. The exponent x is then 
unique, in the sense that it does not depend on the strength of the perturbation H' — H (the coefficient 
A). When this operator is exactly marginal, however, there is a continuum of possible exponents x, which 
depend on the strength of the perturbation. This is exactly what happens in Anderson's original treatment 
of the Fermi gas, and this is why there is a continuum of possible values of the exponent x in that case. 

Both for the "cut-and-glue" and the "swap" quenches, one needs to be able to cut some part of the system 
into two parts. Then the most favorable situation would be a case when the perturbation H' — H ~ A0(.x o ) is 
relevant. Then, a low energy, the degrees of freedom between the part A and B would be effectively decoupled. 
A physical situation where this could be realized corresponds to a Luttinger liquid where one switches on a 
localized potential, which induces some backscattering. The backscattering is known to be always relevant if 
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Figure A3. The two world-sheets that correspond to the "cut-and-glue" quench and to the "swap" quench. 
For the "cut-and-glue" quench, it is a single sheet with a slit. In the case of the "swap" quench, there are 
n sheet corresponding to the n replicas. The top sheet is identified with the bottom one along the dashed 
line. 




(A.18) 
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the interactions between the particles are repulsive. The orthogonality catastrophe exponent x is, however, 
not completely fixed, even in that case. This is because, although the contribution of the backscattering to 
the exponent x is 1/16, as expected in a system with central charge c = 1 (the general result being c/16), 
it is not the only contribution. There is another contribution due to forward scattering , which is always 
marginal, and therefore depends continuously on a parameter (here, the phase shift at the Fermi point). 

On the other hand, one also needs to be able to glue the two parts A and B together, such that the 
total system A U B is homogeneous. In that case, it is the opposite situation which is the most favorable, 
namely the case when <p(x, r) is irrelevant along a line. In a Luttinger liquid, this is the case when the 
interactions are attractive. We thus see that, while the RG flow helps us to realize one of the two operations 
(either to "cut" or to "glue" the two subsystems), it makes the other operation difficult. Again, the "swap" 
quench is even more difficult to implement, as it requires not only to be able to perform the "cut" and "glue" 
operations, but also to apply them simultaneously on the n copies of the system. 



Appendix A. 7. Orthogonality catastrophe from cone and corner free energies on the world-sheet 

The relevance of the world-sheet geometry to the calculation of the entanglement entropy was pointed out 
already in the early work of Holzhcy et al. [10]. Essentially, the Renyi entropy for integer n is nothing but 
the ratio of the partition functions of the field theory on two different world-sheets. One way of seeing that 
is to write the expectation value of the swap operator in the ground state |^>) . . . ® as the overlap 



V/'opi Vrcpi. ) as in Eq. (A. 8), and then use the fact that \ip Iep i.) (resp. Vrcpi )) i s tne ground state of 

Hrepi. (resp. H' rcpl — Sa • H lep i ■ S A ). Since imaginary time evolution with H lcp \ or H^ epl automatically 
acts as a projection onto these states after a long time, one can re-write the overlap as 
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(A.19) 



The states \in) and \out) can be chosen arbitrarily, as long as they have a non-zero overlap with the ground 
state of -ffrepi. and H' lepl respectively. In field theory, the numerator in the right-hand side is precisely 
the partition function on the world-sheet drawn in Fig. A3 (right), while the denominator is the square 
root of the partition function on a world-sheet composed of 2 x n separated strips (note that the square 
root compensates the factor 2). This cartoon is very useful to evaluate the exponent of the orthogonality 
catastrophe associated with the swap, 

|«pl|^-p1.)| ~ (£) ■ (A.20) 

Indeed, the exponent a; swap turns out to be a difference of free energies / = — log Z on the different world- 
sheets. There are extensive parts in these free energies, which scale as L 2 (bulk free energy) and L (surface 
free energy) respectively. However, these contributions appear both in the numerator and the denominator 
of (A.19), and they all cancel. The next leading contribution is of order O(logL) is a conformal field theory; 
it is the free energy / cone that comes from the conical singularity in the world-sheet. The coefficient of 
this term is nothing but the orthogonality catastrophe exponent, as can be seen by taking the logarithm of 
(A.20): 

x logfX/a) ~ -log | (Vfcpi.lVw.)! - f conc - ( A - 21 ) 
It is an interesting exercise in conformal field theory to evaluate / conc for a cone of angle 9; the result was 
derived in a famous paper of Cardy and Peschel [17], and it is given by A (J- — log(L/a). The swap 
with n replicas corresponds to a cone of angle = lixn. This gives the orthogonality catastrophe exponent 
[10, 12] 

' f n--\. (A.22) 



anap 12 V n. 

For the world-sheet corresponding to the "cut-and-glue" quench, the situation is very similar, except that 
there is no cone in the bulk, but rather a corner at the tip of the defect line (Fig. A3). In that case, the 
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orthogonality catastrophe exponent £ cut / g i uo is identified with the coefficient of the O(logL) term in the 
comer free energy, for a corner of internal angle 9 — The evaluation of / corner in a conformal field theory 
was also done by Cardy and Peschel [17]; we recalled this result in section 2.2. The orthogonality catastrophe 
exponent associated to the "cut-and-glue" quench is 

2-cut/gluc — 77T- (A. 23) 



Appendix B. The effect of boundary changing operators 

In this appendix we consider the effect of boundary changing operators on the universal finite-size function 
f(x — i/L), extending the calculation presented in [15] when there are none. We first derive a formula in 
the most general case of four different changes in boundary conditions (Appendix B.l), before specializing 
to more realistic situations (Appendix B.3 and Appendix B.3) relevant to the main text. 



Appendix B.l. Leading universal term 

Let us look at the LBF in imaginary time, when the four possible boundary conditions, denoted by a, /?, 7, 5, 
are different. We have 

T , — of a P~t & f a @ f1 S f aS (TZ -\\ 

■TaPyS — Z JA,B J A ~JB ~JAUB \ D - L ) 

which are represented in figure Bl. 
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Figure Bl. Free energies involved in the calculation of the LBF (Eq. B.l), with four different boundary 
conditions a, /3, 7, 5. Left: f^fg S ■ Middle: f^^B = f'jt +/s*" ^9^ : Iaxjb- The lower boundary is by 
convention at Im(z) = 0. 



In the following we drop all the greek superscripts, to lighten the notations. Let us now consider the 
infinitesimal transformation w w + i5£ inside the strip, and idhjkj otherwise. The corresponding variation 
in the free energy can be expressed using the T yy component of the stress-energy tensor and the boundary 
operators. For /a,b we have 
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and there are similar expressions for /j^^g, Iaub- ^1 ( res P- 4>2, 03, 04) is an operator that changes the 
boundary condition from a to /3 (resp. /3 to 7, 7 to (5, 5 to a). It has dimension hi (resp. /12, ft.3, /14). Notice 
that it is necessary to introduce a cutoff A. Indeed all the free energies defining the LBF are infinite in the 
limit A — > 00, but the crucial point is that the linear combination 
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remains finite. The variation (B.2) can be computed using T yy — T(w) + T (w) and the inverse conformal 
transformation to the upper-half plane, see Fig. B2. The conformal transformation that maps the upper 
half-plane to the asymmetric pants geometry is given by 
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Figure B2. Le/i: Asymmetric pants geometry, with four different boundary conditions a, /?, 7, <5. Right: 
Upper Half plane H, and the four boundary changing operators. 



and the transformation law of the stress tensor reads 

where {w(z),z} is the schwarzian derivative. Using this, Eq. (B.2) becomes 
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The bec operators are inserted at the points Z\ = — 1, Z2 = 0, Z3 = 1/x — 1 and at Z4 = ioo in the z plane. 
Without them u/,(z) = and the calculation boils down to that in Ref. [15]. Vh{z) can be found using the 
conformal Ward identity for chiral operators 
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and the general expression for the chiral four-point function 
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Here z. 



Zj\, H = J2t=i hi and T is a universal function of the anharmonic ratio f. The second 
integral in (B.6) can be deduced from the first by symmetry, replacing x by 1 — x and exchanging the roles 

ai 

c(x 2 - 1)' 



of hi and h 3 . Computing the first integral and repeating the same procedure for f^ B and f^ B cancels 
out all the diverging cutoff-dependent terms, as should be. We get 
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after long algebra. We have introduced the notation 
— > 1 — x ; h\ ^ h 3 
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for the action of replacing x by 1 — x and swapping h\ and /13 in all the preceding terms. After integration 
this gives us the universtal scaling function f(x) for the LBF. The final result reads 
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Of course this function depends on the precise form of T, and f(x) is sensible to the refined details of the 
CFT. A simple case of interest is that of vertex operators 

(V ai (z 1 )V a2 (z 2 )V a3 (z 3 )V ai (z i )}= [] Z T^ ( B - 14 ) 

l<i<j<4 

where the on are the charges, submitted to the neutrality condition JZt—i a i = 0- This form can be identified 
with (B.10), setting hi = af/2 and 

We get 
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The case without any boundary condition changes [15] can be recovered by setting all the charges a; = in 
(B.16), or, equivalently, setting all the dimensions hi — and T(£) = 1 in (B.13). 



Appendix B.2. Subleading semi-universal term 

Using the results presented in Appendix B.l, the calculation of the L _1 logL becomes straightforward. It 
suffice to notice that (T(w)) is modified, in the presence of boundary condition changing operators, to 
(T(ui)0i</>20304) / Then all the arguments given in Sec. 2 apply, and we find 
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which may be evaluated using the Ward identity (B.9). For the pants geometry we find 
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Note that this expression does not depend on the function of the anharmonic ratio T(£). 
Appendix B.3. Special cases 

Special case abba We now specialize to the case abba, with only two boundary changing operators inserted. 
This is relevant to Sec. 3.2.1. We find, setting a 2 = 014 — and a\ = a§ = 2h: 
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Special case abab Here we specialize to the case of 4 bec operators of the vertex type, relevant to Sec. 3.2.2. 
Setting a,\ = a 3 = —a 2 = —u± = y2h, we get 

£r-h) (2X-1+ 2 

.24 J \ 

In the XXZ chain h = h p = =5- (f ) 2 , where R is the compactification radius and S the phase shift at the 
Fermi surface[44]. In the XX limit R = 1 and h p = l/2(p - 1/2) 2 . 
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Appendix C. LBF for an XX chain cut into two equal parts: an exact result 



In this appendix we recall and expand the exact result of [50] , for the LBF of an open XX chain cut in two 
subsystems of the same size. A summary is presented in Appendix C.l. After diagonalization (Appendix 
C.2), we show how the fidelity can be computed as a Cauchy-type determinant, and expressed in closed form 
(Appendix C.4). The asymptotic expansion is explained in Appendix C.5. Finally, we give in Appendix C.6 
two other examples to which the method presented here can be applied. 



Appendix C.l. Summary of the asymptotic expansion 



Before proceeding with the derivation, let us first present the exact result for the asymptotic expansion of 
the LBF. For any integer p > 1, we have 
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The leading terms are given by 
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hp = | (p — 1/2) is the dimension of the (phase shift) boundary changing operator discussed in the text. 
In CFT, 70 is determined from the Cardy-Peschel formula, while 7 i is a consequence of our main result. In 
this particular geometry 7^ is given by 
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for k > 1. We were also able to compute the (non-universal) constant term, which is given by 
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2 Cl 2 (7rp) + (1 - 2p) Cl 2 (2^o) + - Cl 3 (27rp) 



2 log 



G 1 



P 



G 



l + P 



(C.5) 



where C is the Riemann zcta function, G the Barnes G function, and the Cl„ are the n— th Clausen functions 
defined by 



Cl n (x) 



sin(fca;) 



k=l 



cos(kx) 



Ecus^ft 



n even 



odd 



(C.6) 



k=l 



At half filling p = 1/2, the constant can be shown to take the simpler form 

T(3/4) \ 



9 1 
Mo = -x log A + - log 



1 3 C 7C(3) 

— log 7T — + - -Y 



8 TT 



-0.126058180038976 (C.7) 



C is Catalan's constant and A is the Glaishcr-Kinkelin constant. 



Appendix C.2. Diagonalization 

We consider an XX chain of length L, in a transverse magnetic field. The Hamiltonian reads 

Haub = ~J2 (->;., + ^j+i) - h lL ° z i- ( c - 8 ) 
j=i j=i 
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Upon performing a Jordan- Wigner transformation 

< + ia) = c] exp Utt E cjc/j (C.9) 

a] = 2c]c 3 - 1, (CIO) 
it may be rewritten as 

L L 

h A ub = - E {4+^ + c h+i) ~ h J2 { 2c h - *) . ( c - n ) 

and then diagonalized by a standard techniques [38]. Introducing the new set of fermions 

v 1/2 L 

~\ ^sin(j^ m )c], (C.12) 



d U , , , , 

\ /, + 1 / 

with the quasimomenta given by 



E (4+i c i+4 c i+i) E ( 2 ^-l). (C.18) 



771 7T 

<^ = ^7 . m = l,2,...,L, (C.13) 
allows to rewrite Haub as 

L 

^Aub = est + e m d\ n d m , e m = - cos(0 m ) - /i. (C14) 

m— 1 

The ground-state may be obtained by filling the momenta with negative energy: 

\AUB) =d\4...d f N \0) (C.15) 

TV is the number of fermions and depends on h. At zero magnetic field and we have N ~ L/2 (half- filling) . 
In general we assume a constant filling fraction p = N/L. Let us repeat the same procedure for the cut 
Hamiltonian 

H A ® B = H A + H B , (C.16) 

with 

La—1 La 

H * = - E (4+i c i + c h+i) ~ h J2( 2c h - !) ( c - 17 ) 

L L 

#B= - 

j=L A +l ' j=L A +l 

Ha®b is identical to Haub except for the local hoping term — (c^ +\Cl a + ^- c ) that connects subsystem A 
and i? in Haub- The two terms in Eq. (C.16) can be diagonalized separately, introducing the momenta 

®k = 7~~TT ' fc=l,2,...,Ln (C.19) 
for either = ^4 or O = _B. The two new set of fermions are 

^j 1/2 Esin(A A H (C20) 

b l = ir-^)' 2 E si *(A S H (C21) 

The ground state is then given by 

\A®B) = a\a\ . . . a\ A b\b\ ...b\ B |0) (C.22) 

where Nq — pL$i is the number of fermions in subsystem Q. To avoid any unnecessary complications, we 
will always assume pLfi to be an integer. For example at half filling (p = 1/2) this means both La and Lb 
have to be even. 
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Appendix C.3. A simple variation on the Cauchy determinant identity 

We take a small detour to derive a Cauchy-type determinant identity that will be useful to us later. The 
celebrated Cauchy determinant formula is given by 



dct 

l<i ,j <n 



1 



III 



<*<j <n 



<2<J <7l 



(Vi ~ Vj) 



(C.23) 



, x i+VjJ Ul<iJ<n( X i + Vj) 

It can be understood by noticing that both sides of the equality are rational fractions that have the same 
poles and zeroes. In the following we will need to evaluate a determinant with a similar structure. Consider 
a n x n matrix with n even and matrix elements 

(-iy 



Xi + yj 
(-l)j 



1 < i < n/2 



n/2 + l<i<n 



(C.24) 



, x i-n/2 + Vj 

The determinant can be expressed as a product of two Cauchy determinants, using elementary row-column 
manipulations. However, it is much more elegant to evaluate it directly, using a similar reasoning: 

\ 2 r-r(i-j)even 



det {rriij) 

l<i j <n 



2 n / 2 



<i<j<n/2 ( X i 



X j) Ul<i<j<„ iVi - V]) 



(C.25) 



The poles are obvious. Then, the determinant has double zeroes for Xi = xj, i ^ j and simple zeroes for 
Hi = yj, i =/= j provided (— l) 1 = ( — 1)-'. The remaining 2"/ 2 factor is explained by the possible choice of two 
Xi for a given i = 1, . . . , n/2. Despite of its simplicity, we have found no mention of the identity (C.25) in 
the literature. 



Appendix C.4- A closed-form formula for the fidelity 

Using the results of the previous section, the overlap we are looking for can be expressed as a determinant: 

(A ® B\A U B) = (0|aio 2 • • . aw A 6i& 2 • • ■ b NB d\d\ . . . 4,10) (C.26) 

= det (M kl ). (C.27) 

l<k,l<N 

Eq. (C.27) follows from Eq. (C.26) by applying the Wick theorem. The matrix elements are, after some 
algebra, given by 



M kl = (0K4|0) = {-iy z 



skA "k l J l 



M kl = (016^^410) = (-l)'z 



a 



oB 



IB ^k-NA^l 



u 



B 

k~N A 



l<k<N A 



N A <k<N, 



Vl 



with 



[(£ + l)(L n + l)] 



-1/2 



sin6»[ 2 



v2 = cos 62 



(C.28) 
(C.29) 

(C.30) 
(C.31) 
(C.32) 
(C.33) 
(C.34) 



Vl = COS (pi . 

This form allows for fast numerical evaluation using standard linear algebra routines. Even more interesting, 
a remarkable simplification occurs when the two subsystems have the same size La = Lb = L/2: we have 
z A = z B = z, and the same goes for the a^, an d u^, so that we can remove all = A, B superscripts. 
In this case we can use the determinant identity (C.25) of Appendix C.3, and get 



Y\k=Mkf rizii & u.i<k<i<N/2( u k - u if ni<fc<r<wK - «o 



riV 



(fe-/)< 



D 



N 



2 N/2 Z N 



T1 N/2 T1 N , x 



(C.35) 
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Our final result for the LBF therefore reads 



T= -log \(A UB\A®B)\ 2 = -log(D%) 



(C.36) 



Appendix C.5. Extracting the asymptotic expansion 

It is straightforward, if somewhat cumbersome, to extract the asymptotic expansion of Eqs. (C.36) and 
(C.35). Let us explain a possible method on the most complicated term, at general filling fraction p = N/L. 



'n/2 N 

log s= - log i nn^-^) 



fc=i 1=1 

N/2 N 



2 sin 



h + 4>i 



sin 



(C.37) 
(C.38) 



k=l 1=1 

Let us now cut (C.38) in two 

- log B = — log B ice - log B na , (C.39) 

using log(sinx) = log (sine a;) + logx. Since log (sine a; ) is analytic, one can apply twice the Euler-Maclaurin 
formula jj at any order, and check that the expansion takes the regular form 



logS rcg = E 6 2-fc i 



2-fc 



(C.40) 



k=0 



The equality is meant in the sense of asymptotic series. Having performed this expansion we are left with 
the "linearized" term 

N/2 N 

B iiD = n + - <w/ 4 ( c - 41 ) 
A;=i i=i 



2{L + 1) 



N(N+i/2) n/2 t(n +l + 2i- -^j r (N + I - 2i + j^j 



n 



W2)n 



24 
L+2 



(C.42) 



r is Euler's Gamma function. The denominator can be handled at arbitrary order using the Euler-Maclaurin 
formula, and generates only terms of the form (C.40). The prefactor is easy. We are then left with 



N/2 



i=i 

N/2 



with 



C 



= nr(2*-i + / , + 2^)r(iv + i + 2z-^) 

i=l 

= C x D, 



Tr-f G(N + l)G(N + 2)G(^±i)G(^±f±^) 



2-N(N+p/2) G ( N + 2 ) G ( A^3) G (£±l) G (£±2) 



"/* r (2* - 1 + p + 2^^) r (n + i + 2% 



2i 
L+2 



i = l 



T (2i - 1 + p) T (N + 1 + 2i) 



(C.43) 

(C.44) 
(C.45) 

(C.46) 
(C.47) 



tt El„ f(i) ~ fa /(*) dx + s{a) + s{b) + EfcLi Iff [/ (2fc_1) (6) - /t 2 *" 1 ' (a)] where /( n ) is the n-th derivative of / and the 
B2k are the Bernoulli numbers, defined through the exponential generating function x(e x — l)~ 1 = Efe^Li ^W xk - The nrs * ^ ew 
are given by B 2 = 1/6, B 4 = -1/30, B 6 = 1/42. 
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where we have introduced the Barnes G-function. G can be obtained by induction, using the functional 
relation G(z + 1) = T(z)G(z) and the duplication formula T(z)T(z + 1/2) = 2 1 - 2z % / 7 rr(2z). This is useful 
because the asymptotic expansion of G{z) for large z is known: 



logG(z + 1) = \z 2 log* - \z 2 + l ^±z - i logz + 1 - log^ + f) 2* , (C.48) 

fe=i 

where A is the Glaisher-Kinkelin constant. The ^fe are known but unimportant to us. Hence 

oo 

- log C = (c 2 L 2 + Cl L + c ) log L + J2 c 2- k L 2 ' k (C.49) 

fc=0 

The last remaining term is 

-log£» = -Vlog l-± T L+2 J . V , ^ (C.50) 

8 ^ I r (2 2 - 1 + P ) r (iv + 1 + 2i) J 1 ; 

We wish to use once again the Euler-Maclaurin formula. Because of the L + 2 denominators, it is much more 
convenient to work as a function of M = L + 2. It is also easy to see that the integration will generate a 
denominator proportional to M — 1. The integral of logr(u) is the generalized polygamma function ip^~ 2 '(u). 
Its asymptotic expansion is simply obtained by integrating the Stirling formula 

log T(z + 1) = z log z - z + - log(2,z) + J2 (C51) 

We get 

oo 

- (M - 1) logD = (d 2 M 2 + dxM + d + d^M' 1 ) logM + ^ d' 2 _ k M 2 ~ k . (C.52) 

fc=0 

Notice that the power series prefactor of logAf terminates at M _1 . However, it does not in terms of L. 
Indeed expanding M^ 1 = {L + 2)" 1 and (M - l)" 1 = (L + in power series, we get 

/ oo ,_ , k \ oo 

- logD = rf 2 L + 3d 2 + di + K 2 * -1 ~ X ) d -i - do - di - rf 2] logL + d2-kL 2 ~ k , (C.53) 

V fe=l / fc=0 

with, in particular, d 2 = p 2 /2, d t = -3/9 2 /2, d = (l-6p+lV)/8 and d_ x = (-l + 6p- 6p 2 )/6. Eq. (C.53) 
already gives all the L - ' logL (fc > 1) terms advertized in Appendix C.l. In principle it is also possible to 
keep track of the L 2 ~ k prefactors, although these become rapidly gruesome. We were able to do so up to 
order k = 2. In the end — In B admits the asymptotic series 

/ oo \ oo 

- Infl = [J2 bi-kL 1 -" log L + ]T b' 2 _ k L 2 - k (C.54) 



with 



\k=0 / fc=o 



b x =0 (C.55) 

b a = (1 + 6p - 6p 2 ) /24 (C.56) 

fe_ p = (-l) p+1 (2 p+1 -l)(l-6p + 6p 2 )/24 , P >1 (C.57) 
W 1 2 i 9 i C(3)-Cl 3 (27rp) 

fo 2 =2^ l0 S 2 + ^2 ( C - 58 ) 



b , = 3(C(3)-C1 3 (2^)) | (2-3p)Cl 2 (2irp)-3Cl 2 M 

1 27T 2 27T 



(3-4p)Cl 2 (2^ j o) -4Cl 2 (7rp) 24C(3) - Cl 3 (27rp) 23 - 78p + 54p 2 . 
fo o = ^ + ~ 2 S logsm.p 

(l-2,)logsinf - 1 + 3 + 2 7 ^ log2 + 5 - 18 2 f ^ log, + log [A 2 G (l + f ) G ■ (C.60) 
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Performing the same expansion on all the other terms in (C.35), we observe that all the terms proportional 
to L 2 , LlogL, L cancel each other, as should be. We finally obtain 



fe=0 



with the 7fc given by 
7o = 



7/c log L + [i k 



7fe = 



( 2 *+i - 1) /i 



- 3/i, 



k > 1 



(C.61) 

(C.62) 
(C.63) 



and no by Eq. (C.5). Notice that it is possible to put the L fc logL corrections in a more compact form. 
Introducing the effective length 



L K = 



(L + 1)(L + 2) 



£ + 2/3 



the series expansion reads 



T{L e ) = 70 log L e + fi Q + 7i 



log£ e 



E 

fe=i 



(C.64) 



(C.65) 



Appendix C.6. Related examples 

It is possible to derive similar exact formulae in two other geometries, that we present here. 



Open-periodic overlap: The first possibility is to consider the overlap between two full chains of length L, 
one with open boundary conditions, and the other with periodic boundary conditions. We have 



(open|periodic) 



det . , 

i<k,i<N I y L(L + 1) 



lir ' 


cos 


'(2fc-l-AT)7T mn ~ 
2L 1 




L + l_ 


COS 


'(2fc-l-jV)7r" 
L 


— COS 


L+l 



(C.66) 



This determinant can also be put in the form (C.24), using the symmetry k — > N + 1 — k. Evaluating this 
in explicit form, and performing an asymptotic expansion similar to that described in Appendix C.5, we get 



log |(open|periodic)|^ 



1 



h p log L + no 



1 



log L ^ // fe ' 



L 



E 

fe=i 



1 \ °° 



k=2 



; logi 



(C.67) 



Fidelity in a periodic system: This was the other situation considered in [15]. When the periodic system of 
length L is cut into two open chains of size L/2, it is also possible to put the fidelity under the form (C.24). 
After asymptotic expansion we find 



P(L) 



2hp)\ogL-Ah p L -r { . 



^E(- 2 )^T + E^ (C-68) 

fe=2 k=0 

Notice that the L~ x \ogh correction vanishes at half-filling, contrary to the other two examples. 
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